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SUMMARY 


This  report  describes  the  work  performed  under  Contract 
AF  04(611) -9073  during  Yhe  third  quarter  of  the  contract  period, 
beginning  December  1963>^The  objectives  of  the  study; mm* 

(a)  To  develop  design  tools,  in  the  form  of  digital  computer 
programs,  for  the  nozzle  designer  to  use  for  calculating  the 
behavior  of  materials  exposed  to  a  solid-propellant  rocket-exhaust 
environment. 

(b)  To  marshall  existing  theoretical  methods,  and  tecr-ssALgmj 
or  develop  new  methods  wi»®*C3*58JBSSSPp,  to  gain  a  better 


understanding  of  the  fundamental  mechanisms#  a  a  soei  atgfT"vy;ffcrt  cxfroggcT 
material*  bilia v  i uirDL:a,' 

-  - - ~~  ujW-7 

#ork  during  rf  8B3E8BSE^BasA  directed  toward 

development  of  two  digital- computer  programs.  O ne-»^««w»( 


applicable  to  materials  such  as  graphite  and  tungsten,  that, 
gemwBWf  erode  with  one  moving  boundary  or  do  not  erode  at  all. 

The  other  applicable  to  the  reinforced  plastic  materials  that 
form  a  char  layer.  effort  in  the 

areas  of  transient  conduction,  equilibrium  chemistry,  reaction 
kinetics,  and  mass  transport,  as  applied  to  the  solid-propellant 
rocket  problem,  is  described  in  detail. 

Further,  an  arc-plasma  generator  rocket-simulation  test  pro¬ 
gram  has  been  conducted  to  study  the  effects  of  chemical  and 
mechanical  erosion  of  several  material  types,  and  the  results  of 
43  nozzle  firings  are  reported.  A  second  test  program  to  investi¬ 
gate  the  erosive  effects  of  Al^O^  particle  impaction  has  been 
formulated,  and  the  pror osed  test  plan  is  described.  ^  j  blit-  flu- 
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A  STUDY  OF  SOLID-PROPELLANT  ROCKET 
MOTOR  EXPOSED  MATERIALS  BEHAVIOR 

1.  INTRODUCTION 

During  the  past  several  years,  significant  theoretical  tech¬ 
niques  have  been  developed  by  various  investigators  for  character¬ 
izing  many  of  the  complex .  individual  phenomena  occurring  in  the 
ablation  process  of  solid-prcpeiianu  rocket  wall  materials.  These 
techniques  have  been  used  to  gain  a  better  understanding  of  abla¬ 
tion  phenomena  under  simplified  conditions  and  to  guide,  to  a 
certain  extent,  the  design  of  aft-closures  and  nozzles.  The  proven 
utility  of  the  theoretical  approach  in  this  area  has  made  meaning¬ 
ful  the  possibility  of  developing  an  integrated,  more  general,  ana¬ 
lytical  technique  for  the  design  of  solid-propellant  rocket  wall 
materials  and  configurations.  The  need  for  such  a  technique  is 
becoming  more  urgent  as  motors  increase  in  size,  and  thus  render 
the  traditional  cut-and-try  approach  more  costly  and  time-consuming. 

Recognizing  this  need,  the  objectives  of  this  study  are  two¬ 
fold: 

(1)  To  develop  design  tools  for  the  nozzle  designer  to  use 
for  calculating  the  behavior  of  materials  exposed  to  a  solid- 
propellant  rocket -exhaust  environment.  These  design  tools  will 
be  in  the  form  of  digital  computer  programs. 

(2)  To  marshall  existing  theoretical  methods,  and  to  extend 
these  or  develop  new  methods  where  necessary,  to  gain  a  better 
understanding  of  the  fundamental  mechanisms  associated  with  exposed 
materials  behavior  in  a  rocket-nozzle  environment.  This  second 
objective  is,  of  course,  coupled  with  the  first  in  that  an  appre¬ 
ciation  of  the  fundamental  ablation  mechanisms  is  required  to  calcu¬ 
late  their  effects. 

Both  theoretical  and  experimental  studies  are  being  conducted 
to  accomplish  the  objectives  of  the  program.  The  relations  charac¬ 
terizing  the  ablation  process  are  being  formulated  largely  from 
theoretical  considerations,  but  also  draw  on  experimental  work  in 
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those  areas  not  amenable  to  theoretical  treatment.  Existing  test 
data  are  being  used  where  available,  and  small-scale  laboratory 
tests  are  being  conducted  as  a  part  of  this  program  to  provide 
information  in  several  areas  where  it  is  required.  The  flow  chart 
of  Figure  1  indicates  how  the  various  areas  of  investigations  are 
interrelated. 

As  indicated  in  Figure  1,  the  theoretical  studies  are  divided 
into  two  major  categories,  the  first  having  to  do  with  materials 
that  erode  with  one  "moving  boundary"  such  as  graphite  (nonablating 
refractories  are  a  special  case  of  this  category) ,  and  the  second 
with  materials  such  as  the  reinforced  plastics,  or  composites, 
whose  ablation  can  be  characterized  by  two  moving  boundaries.  The 
experimental  studies  are  being  carried  out  in  the  Vidya  1-megawatt 
arc-plasma  generator  facility;  they  consist  of  graphite  and  composite- 
nozzle  ablation  tests  and  metallic-oxide  particle-impaction  tests. 

During  the  third  quarter  of  the  program,  the  theoretical  effort 
continued  in  the  following  areas:  thermal  behavior  of  nonmelting, 
noncharring  materials,  thermal  behavior  of  charring  materials,  and 
wall  chemical  erosion  (both  diffusion  and  kinetically  controlled) . 

The  accomplishments  here  are  reported  in  Section  2.  The  experi¬ 
mental  effort  durign  the  quarter  consisted  of  conducting  and 
interpreting  43  nozzle  firings  under  the  chemical-erosion  program 
and  formulating  test  plans  for  the  particle  impact  program.  This 
work  is  discussed  in  Section  3. 

2.  THEORETICAL  STUDIES 
2.1  Introduction 

For  the  purpose  of  classification,  all  engineering  effort  on 
the  program  that  does  not  directly  involve  an  experimental  test 
program  is  called  "theoretical.  11  This  is  actually  a  misnomer, 
since  the  effort  is  based  largely  on  proven  physical  and  mathe¬ 
matical  relations,  and  is  of  a  very  practical  and  applied  nature; 
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that  is,  the  effort  is  directed  toward  the  construction  of  calcu¬ 
lation  tools  for  characterizing  the  phenomena  occurring  at  an 
eroding  nozzle  wall.  Bees  jse  many  physical  relations  are  needed 
for  this  characterization,  and  because  the  interactions  between 
them  are  quite  complex,  an  analysis  flow  diagram  was  prepared 
(Pig.  2)  to  help  clarify  the  role  of  each  area  of  study. 

Figure  2  is  itself  necessarily  complex,  but  it  is  hoped  that 
it  will  be  helpful  in  placing  each  phase  of  the  effort  discussed 
herein  in  proper  perspective.  The  row  of  boxes  at  the  top  of  the 
figure  represents  the  various  inputs  required  to  perform  nozzle 
wall  erosion  calculations;  except  for  the  erosion  effects  of 
particle  impact,  none  of  these  specific  data  are  being  sought 
experimentally  under  the  present  contract.  The  program  emphasis, 
rather,  is  on  the  "calculation"  area  of  the  figure,  to  provide  a 
technique  for  determining  nozzle  erosion  rates  and  wall  tempera¬ 
ture  from  the  input  data. 

The  work  reported  in  this  section  is  in  the  areas  of  surface 
chemical  reactions  and  internal  material  behavior.  As  che  reader 
reviews  the  work,  it  is  recommended  that  he  refer  to  Figure  2 
occasionally  to  review  the  role  of  each  study. 

2.2  Thermal  Response  of  Materials  Which  Ablate  With  One  Moving 
Boundary 

2.2.1  Introduction 

The  class  of  materials  considered  here  is  that  class  which 
undergoes  physical  and  chemical  change  only  at  the  surface  which 
is  exposed  to  the  exhaust  environment.  The  most  common  example 
of  these  materials  encountered  in  rocket  motors  is  graphite.  A 
prediction  of  the  response  of  this  class  of  materials  to  an  exhaust 
environment  involves  a  consideration  of  the  internal  thermal 
response  (through  conduction  of  heat  away  from  the  exposed  surface) 
and  the  thermal  and  chemical  interaction  between  the  exposed  sur¬ 
face  and  the  exhaust  stream. 
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2.2.2  Surface  chemical  reactions 

A  study  of  the  chemical  reactions  occurring  at  the  exposed 
surface  has  a  two-fold  purpose,  first,  to  determine  the  rate  at 
which  wall  material  is  consumed  by  chemical  reactions  and,  second, 
to  enable  specification  of  the  energy  state  (enthalpy)  of  the 
mixture  of  gases  existing  at  the  wall.  In  general,  these  quantities 
are  determined  as  functions  of  temperature  and  exhaust  stream  param¬ 
eters  (pressure,  transfer  coefficients,  etc.)  and  are  needed  to 
supply  a  boundary  condition  to  the  calculation  of  the  temperatures 
in  the  interior  of  the  wall  material. 

The  basic  analytical  technique  used  to  calculate  the  desired 
quantities  is  a  surface  mass  balance.  This  mass  balance  simply 
states  that  the  net  rate  of  arrival  of  an  identifiable  molecular 
species  to  the  surface  must  vanish;  that  is,  the  sum  of  the  rates 
of  supply  of  a  molecular  species  due  to  gas  convection  normal  to 
the  surface,  diffusion  through  the  boundary  layer,  condensed  phase 
mass  addition  and  chemical  reactions  at  the  surface  must  be  zero. 
Further  discussion  of  the  mass  balance  technique  can  be  found  in 
References  1  and  2. 

Consider  the  control  surface  shown  in  the  following  sketch; 
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(pv)w 


* 


roa  Kja  +  *d  Kjd  +  wj  (1) 


where 

(pv)w  net  convective  mass-velocity  normal  to  the  surface 


mass  fraction  of  species  j  in  the  surface  gas 
mixture 


D.  diffusion  coefficient  of  species  j  through  the 

mixture 

.* 

n»r  mass  rate  of  removal  of  condensed  phase  material  from 

A  the  surface 

* 

K.  mass  fraction  of  species  j  in  the  condensed  phase 

-*  material  being  removed 

m,  mass  rate  of  erosion  of  virgin  ablator 


mass  fraction  of  species  j  in  the  ablation  material 

mass  rate  of  deposition  of  condensed  phase  material 

mass  fraction  of  species  j  in  the  condensed  phase 
material  being  deposited 


net  rate  of  production  of  species  j  due  to  surface 
chemical  reactions 


The  chemical  production  term  can,  in  general,  be  expressed  in 

the  form  of  a  reaction-rate  term  multiplied  by  a  difference  in 
partial  pressures.  For  example,  if  the  reaction 

C*  +  C02  £  2C0 


is  the  only  reaction  involving  C02  at  the  surface,  then 


(2) 


*  =  k  ( 

C02  f  Up 


where 

Jcf  reaction  rate  constant 
Pcq  partial  pressure  of  CO 
pCO  partial  pressure  of  CO. 


Kp  equilibrium  constant  for  reaction  (2) 

The  reaction  rate  kf  is  a  function  of  temperature,  as  is  K  , 
and  the  Pj  are  related  to  the  through  the  equation 

Mipi 
K  =  — 1—1 

J  PM 

where  p  is  the  system  pressure  and  M  the  molecular  weight  of  the 
gas  mixture  at  the  surface.  Writing  Equation  (1)  for  each  species 
considered  yields  J  equations  which  are,  in  principle,  indepen¬ 
dent.  It  will  be  shown  later  that  in  a  practical  sense  they  may 
not  all  be  independent.  We  have,  in  addition,  the  equations 


I  V 


-  Z- 


j  P 


Httnce,  we  have  J+2  equations.  However, 


the  unknown  quantities  are 
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Ky  dKj/by,  (pv)^,  and  M,  so  that  there  are  2J+2  unknowns.  We 
might  introduce  the  boundary- layer  equations,  but  this  introduces 
considerable  complexity  quickly.  Instead,  we  would  like  to  intro¬ 
duce  a  general  transport  equation,  which  really  represents  an  inte¬ 
grated  form  of  the  boundary-layer  equations,  in  the  form 

-<^)w  -  PeueCMj<Kjw  '  V  (6) 

Equation  (6)  can  only  reasonably  be  expected  to  hold  for  a  non¬ 
reacting  boundary  layer.  If  reactions  occur  in  the  boundary  layer, 
its  application  to  a  reacting  species  can  be  shown  to  result  in  a 
logical  contradiction,  at  least  under  some,  if  not  under  ail  circum¬ 
stances.  The  idea  of  introducing  Equation  (6)  is  still  very  appeal¬ 
ing,  however,  because  if  the  K.  are  known,  the  number  of  equations 

J  ® 

in  the  system  is  increased  to  2J+2  without  increasing  the  number 
of  unknowns,  and  we  have  immediately  a  determinate  system.  It  will 
be  shown  in  the  next  two  sections  how  Equation  (6)  can  be  applied 
to  many  situations  of  interest,  even  without  the  restrictive  assump¬ 
tion  of  a  chemically  frozen  boundary  layer.  Because  some  simplifi¬ 
cations  result  when  the  gas  mixture  at  the  surface  is  in  equilibrium, 
we  shall  consider  that  case  first  before  turning  to  the  general 
situation  with  finite  surface-reaction  rates. 

2. 2. 2.1  Surface  equilibrium 

Equilibrium  at  the  surface  will  be  achieved  when  the 
reaction  rates  are  very  large.  Consider  the  general  chemical 
reaction 


rNl  +  nNj  £  qNm  +  sNk  (7) 

where  is  a  general  chemical  symbol  and  r,  n,  q,  and  s  are 
the  stoichiometric  coefficients.  Ideally,  the  rate  of  production 
of  species  j  by  this  general  reaction  can  be  expressed  as 
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w . 
J 


s 

pk 


PJ 


(8) 


where  K  j  is  the  equilibrium  constant  for  the  general  reaction. 
This  will  be  taken  to  represent  both  surface  and  gas-phase  reactions. 
If  the  reaction  given  by  Equation  (7)  is  nonideal,  the  exact  form 
of  Equation  (8)  may  be  altered.  For  example,  the  exponents  on  the 
partial  pressures  may  not  equal  the  stoi  .iometric  coefficients. 
However,  the  net  production  of  any  species  is  still  expressed  as 
the  difference  between  a  forward  and  a  reverse  reaction  rate. 

If  we  substitute  Equation  (8)  for  in  Equation  (1)  and 

divide  by  kf  ^  we  have 


Now  we  notice  that  the  J  unknowns,  dK^/dy,  do  not  appear  and, 
further,  Equation  (10)  is  exactly  the  equilibrium  relation  among 
the  partial  pressures  of  the  species  involved  in  reaction  (7) .  It 
should  be  noted  that  the  ideal  form  chosen  for  the  general  reaction 
rate,  Equation  (8) ,  in  no  way  influences  the  conclusion.  The  same 


-9- 


conclueion  is  reached  if  one  recognizes  that  only  at  equilibrium 
the  forward  rate  equals  the  reverse  rate. 

Por  those  molecules  which  do  not  react  with  the  surface 
but  which  may  experience  gas  phase  reactions,  the  same  results  can 
be  obtained  by  writing  a  mass  balance  for  a  small  but  still  finite 
volume  element  adjacent  to  the  surface.  As  k^  -*■  0,  all  of  the 

terms  in  the  mass  balance  will  vanish  except  for  the  equilibrium 
relation.  Equation  (10). 

We  have  not,  however,  completely  abrogated  the  difficulty 

with  the  system  of  equations  resulting  from  Equation  (1).  We  have, 

to  be  sure,  reduced  the  number  of  unknowns  to  J+2 ,  but  it  can  be 

shown  that  of  the  equilibrium  relations  among  J  molecular  species 

only  J-I  are  independent,  where  I  is  the  number  of  elements 

present  in  the  system;  thus,  we  have  only  J+2  -  I  independent 

equations,  and  J+2  unknowns. 

To  supply  the  additional  equations  required  we  shall  make 

a  mass  balance  on  the  chemical  elements  at  the  surface.  If  j 

molecules  involve  element  i  and  a. ,  is  the  mass  fraction  of 

*  J 

element  i  in  molecule  j ,  we  can  multiply  Equation  ( 1  )  by  ai j 
and  sum  over  j,  obtaining 


j 


where  the  tilde  indicates  an  elemental  mass-fraction.  Since  elements 
are  conserved. 


I  aij  *j E  0 
j 
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If  the  are  equal  or  if  we  accept  a  weighted  average  diffusion 

coefficient.  Equation  (11)  becomes 


<pv)w 


=  ft 


a 


K. 

la 


+  *.* 


(12) 


Lees,  Reference  3,  presents  a  plausibility  argument  suggesting  that 
the  diffusional  mass  flux  of  the  elements,  which  are  of  course  con¬ 
served  through  the  boundary  layer,  can  be  expressed  as 


o  u  CM 
e  e  M. 

i 


(K. 

ie 


-  K.  ) 
iw 


(13) 


Then,  Equation  (12)  becomes 


(pv,w  +  PeVV  <*iw  -  +  K  *ir  ■  *a  «ia  +  *id  <14> 


If  values  are  available  for  p  u  C  ,  A*,  and  K.  ,  the  elemental 

e  e  d  ie 

composition  of  the  gas  mixture  at  the  wall  can  be  determined  as  a 

function  of  wall  erosion  rate,  ft  .  For  example,  if  ft*  =  ft*  =  o, 

a  r  d 

Equation  (14)  assumes  the  familiar  form  (see  Ref.  1) 


K. 

iw 


K.  +  B'K. 
ie _ ia 

1  +  B' 


(15) 


where 


B* 


(pv)w 


It  is  not  necessary  for  the  Cu  to  be  equal  in  order  to  get  a 

Mjl 

solution,  but  it  certainly  simplifies  computation. 

Then  the  equilibrium  chemical  system  previously  developed 
becomes  determinate  with  the  addition  of  I  equations  of  the  form 
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Y  a.  .  K,  =  K. 

L>  ij  jw  iw 
j 


(16) 


2, 2, 2. 2  Pinite  surface  reaction  rates 

If  some  of  the  surface  reactions  are  slow  enough  that  the 
terms  on  the  right-hand  side  of  Equation  (9)  are  significant,  we 
are  still  caught  in  the  original  equation  shortage  dilemma.  The 
escape  here  is  to  devise  a  model  for  which  the  transport  equation. 
Equation  (6),  might  reasonably  be  expected  to  apply,  and  then  to  use 
it  to  supply  the  additional  equations  needed  to  arrive  at  a  determi¬ 
nate  system. 

We  shall  consider  a  subsystem  consisting  of  a  group  of 
molec  lar  species  which  share  a  common  element.  The  common ,  or 
bound,  element  is  distributed  among  the  species  of  the  subsystem 
according  to  the  dictates  of  overall  system  equilibrium,  but  enters 
or  leaves  the  subsystem  only  bv  wav  of  specified  surface  reactions. 
The  rates  of  the  surface  reactions  involving  the  bound  element  can 
be  anything  from  zero  to  infinity.  All  the  other  elements  in  the 
subsystem  may  enter  or  leave  only  by  way  of  fast  reactions. 

A  mass  balance  on  the  bound  element  at  the  surface  is 
(assuming  no  condensed  phase  transport) 

‘pv'w  *?•  -  dL  (w^)v  "  *a  +  *i  (17) 

where  the  superscript  b  refers  to  the  bound  element  alone  in  order 
to  distinguish  these  atoms  from  other  atoms  of  the  same  element 
which  may  be  present  outside  the  subsystem. 

Since  the  bound  element  is  conserved  in  the  boundary  layer, 
we  might  expect  that  the  transport  equation.  Equation  (6),  would  be 
valid.  Then  we  have 


.  / 1 1* 
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P  u  CM 
re  e  M. 
1 


£.b) 

le 


(18) 


If  we  associate  a  subsystem  with  each  of  the  "slow"  surface  reactions, 
the  bound  element  balance.  Equation  (17),  replaces  the  surface 
balance  on  the  reacting  molecular  species,  Equation  (1) .  Hence, 
the  concentration  gradient  of  the  reacting  species  does  not  appear 
as  an  unknown  in  the  system  of  equations,  but  we  have  introduced 
K.  and  (dK./dy)  as  additional  unknowns.  One  of  the  additional 

XW  1  W 

equations  needed  is  supplied  by  Equation  (18) ;  the  other  comes  from 
conservation  of  the  bound  element  and  is 


1 


a  b  K. 
ij  J 


(19) 


where  the  summation  is  only  over  those  molecules  in  the  subsystem. 
This  concept  of  the  equilibrium  microcosm  containing  a  conserved 
element  was  conceived  by  Kendall  and  first  presented  in  a  paper  by 
Kendall  and  Rindal,  Reference  4.  The  principal  utility  of  the  con¬ 
cept  is  to  allow  the  use  of  the  macroscopic  equations  to  quantify 
mass  diffusion  to  a  general  reacting  surface.  Inherent  in  the 
definition  of  the  subsystem  is  that  gas-phase  reactions  affecting 
the  species  in  the  subsystem  are  fast.  There  appears  to  be  no  great 
difficulty  in  relaxing  this  assumption  and  this  will  be  pursued  in 
future  effort. 

To  illustrate  this  concept,  consider  a  subsystem  composed 
of  acetylene  and  methane.  The  bound  element  is  taken  to  be  carbon 
and  enters  the  subsystem  only  through  the  reaction 


2C*  +  Hj 


C2  H2 


(20) 


Then,  since  the  mass  fraction  of  carbon  in  C2H2  is  24/26>  the 
production  rate  of  bound  carbon  is 


-13- 


(21) 


where  and  are  the  reaction  rate  constant  and  equilibrium 
constant,  respectively,  for  reaction  (20).  The  statement  that  the 
bound  element  distributes  itself  in  the  subsystem  according  to  the 
dictates  of  equilibrium  means  that 


CPcO* 

X  f  Pit 

C?2  V^H2  J 


I 

where  K 

P 


is  the  equilibrium  constant  for  the  reaction 


C2H2  +  3H2  £  2CH4 


The  mass  balance  on  bound  carbon  at  the  surface  becomes 


+  Peu«CH  <*c£  '  *Ce> 


24  /  H2 

26  kf  (  K  ’  PC  H9 

P  4.  Z 


(22) 


(23) 


And  the  conservation  equation  is 


b 


24  12 

26  C«H0  16  CH. 

2  2  4 


(24) 


The  specification  of  the  state  of  the  subsystem  then  requires  the 
simultaneous  solution  of  Equations  (22),  (23),  and  (24). 


1 
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2. 2. 2. 3  Graphite 'Surface  reactions 

For  the  graphite  wall  in  an  H,  C,  N,  0,  Cl  exhaust,  three 
such  subsystems  were  chosen.  They  were 

(a)  The  compounds  C2H2,  CH^,  C2H,  C^H,  C^H,  C^-H,  CH2,  C^H2, 

and  C  H  ,  containing  bound  carbon  produced  by  reaction 
(20)  ; 

(b)  The  compounds  H20,  OH,  O,  and  02,  containing  bound 

oxygen  produced  (in  this  case  consumed  is  a  better 
word)  by  the  reaction 

C*  +  H2Ob  £  CO  +  H2  (25) 

(c)  The  compound  C02  alone  produced  by  the  reaction 

C*  +  C02  £  2C0  (26) 

The  remaining  compounds  considered  were  assumed  to  satisfy  the  con¬ 
ventional  equilibrium  relations  among  their  partial  pressures. 

„It  should  be  apparent  that  there  is  some  latitude  in  the 
selection  of  the  molecules  which  are  included  in  the  various  sub¬ 
systems.  For  example,  C  and  02  could  well  be  included  with  the 
CO,,  in  the  third  subsystem,  rather  than  in  the  second  subsystem. 

In  principle,  one  can  also  consider  two  (or  more)  parallel  reactions 
as  contributing  to  the  supply  (or  consumption)  of  the  "bound" 
elements  in  the  subsystem.  The  choice  of  the  three  specific  reactions 
for  the  graphite  wall  was  prompted  by  the  availability  of  data  and 
the  significance  of  molecular  concentrations.  For  example,  data  on 
the  reaction 

2C*  +  2 OH  ^  2C0  +  H2 

were  lacking,  and  very  little  free  02  is  present  in  the  exhaust  so 
that  the  reaction 
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2C  +  02  2C0 

was  not  considered.  The  final  test  of  the  adequacy  of  the  choice 
of  member  compounds  in  the  various  subsystems  will  be  in  how  well 
the  calculation  fits  observed  performance.  All  of  the  molecular 
species  which  react  with  the  wall  to  any  extent  should  in  some 
manner  be  tied  to  a  reaction  which  is  kinetically  controlled,  since 
this  will  insure  compliance  with  the  observed  phenomenon  of  zero 
erosion  at  low  wall  temperatures.  Thus,  even  though  the  reaction  of 
OH  with  graphite  was  not  considered  directly,  the  inclusion  of  OH  in 
subsystem  (b)  ties  that  molecule  to  a  kinetically  controlled  wall 
reaction. 

For  acetylene  production,  the  reaction  rate  was  determined 
by  an  extrapolation  of  data  given  by  the  Union  Carbide  Research 
Institute,  References  5  and  6.  Expressed  in  an  Arrhenius  form,  the 
rate  is  given  by 

OCf)c  H  =  4- 525x10  exp  (=&■ j,2-0-0^)  — - ~ -  (27) 

22  V  y  ft  -sec-atmosphere 

where  T  is  in  °K. 

For  the  CO2  ~  graphite  reaction  given  by  Equation  ( 2) ,  the 

rate  was  taken  to  be 

(kf) co  =  13. 83xl07  exp  C'-2-^0^)  - - — - -  (28) 

2  V.  y  ft  -sec-atmosphere 

where  T  is  in  °K.  This  represents  an  extrapolation  of  an  equation 
given  by  Walker,  Rusinko,  and  Austin,  Reference  7,  which  w»s  based 
on  measurements  by  Gulbranson  and  Andrew,  Reference  8,  and  Armington, 

Reference  9. 

Blyholder  and  Eyring,  Reference  10,  es-oress  the  rate  of 
reaction  (25)  as 
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where  K  is  the  equilibrium  constant  for  reaction  (25) .  The 

Jr 

rate  constants  are  given  by 


— - -  (30) 

ft  -sec-atmosphere 

and 

=  4 024^  4l9,,830N\  - lb -  (31) 

Ko  T  J  exp  T  J  atmosphere 

where,  again,  T  is  in  °K. 

Using  the  rate  constants  given  by  Equations  (27),  (28),  (30), 
and  (31),  the  nondimensional  erosion  rate  of  a  graphite  wall  in  a 
typical  exhaust  was  calculated.  Figure  3  is  a  plot  of 


t 

B 


(PV)W 

p  u  C„ 
re  e  M 


as  a  function  of  wall  temperature.  For  comparison  purposes.  Fig¬ 
ure  3  also  shows  the  predicted  erosion  if  all  surface  chemical 
reactions  are  very  fast,  that  is,  if  equilibrium  is  achieved.  The 
profound  effect  of  reaction  kinetics  is  clear.  On  the  basis  of 
very  sparse  nozzle  erosion  data,  it  appears  that  the  kinetically 
controlled  prediction  is  closer  to  the  real  situation  than  is  the 
prediction  based  on  fast  reactions. 

In  order  to  predict  with  accuracy  the  performance  of  a 
graphite  wall  in  a  rocket  exhaust,  one  must  have  reliable  rate  data 
for  the  particular  type  of  graphite  and  graphite  reactions  involved. 
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2.2.3  One-rooving-boundary  computer  program 

A  digital  computer  program  has  been  developed  to  calculate 
the  thermal  response  of  materials  which  erode  with  one  moving 
boundary.  There  are  two  fairly  distinct,  but  coupled,  parts  to  the 
program.  One  deals  with  the  conduction  of  heat  to  the  interior, 
the  other  with  the  convection  of  heat  and  mass  to  the  ablating  sur¬ 
face. 

The  internal  thermal  response  portion  of  the  program  is  an 
explicit  finite  difference  solution  of  the  differential  equation  of 
heat  conduction.  The  form  and  derivation  of  this  equation  are 
discussed  at  length  in  References  1  and  7  and  will  not  be  repeated 
here.  The  differential  equation  is  transformed  into  a  coordinate 
system  fixed  to  the  ablating  surface  as  follows.  The  basic  body 
shape  considered  is  two-dimensional  and  axi symmetric,  that  is,  it 
can  be  specified  in  r-  and  z-coordinates.  The  differential  equa¬ 
tion  for  the  conduction  of  heat  in  axisymmetric  cylindrical  coordi¬ 
nates  is  transformed  from  r-,  z-,  and  ©-coordinates  to  x-;  z’-,  and 
6 1 -coordinates  via  the  transformations 

x  =  r  -  r 

8 

z  *  -  z 

0'  =  0 

Where  r#  is  the  local  surface  radius.  Because  the  surface  can 
recede,  r  is  a  function  of  both  axial  position,  z,  and  time,  0. 

o 

The  following  sketch  illustrates  the  basic  body  shape  and 
nodal  positions: 


Grid  lines  of  the  coordinate  system  are  constant-x  and  constant-z 
lines.  Nodes  are  positioned  at  the  intersections  of  the  grid  lines. 
The  position  of  any  node  is  specified  by  the  double  index  m,  n. 

The  axial,  or  z-position  index,  m,  runs  from  1  to  ^max*  The  radial, 
or  x-position  index,  n,  runs  from  1  to  nmax  and,  in  general,  nfflax 
is  different  for  each  value  of  m. 

The  inner  surface  shape  is  arbitrary  so  long  as  it  is  con¬ 
tinuous.  The  selection  of  Az  and  Ax  is  arbitrary  within  the 
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context  of  the  stability  of  the  explicit  solution  to  the  difference 
equations . 

Of  course,  the  spacing.  Ax,  between  adjacent  rows  of  nodes 
is  constant  along  any  row,  and  the  spacing,  Az,  between  adjacent 
columns  of  nodes  is  constant  up  any  column,  but  either  or  both  may 
vary  from  row  to  row  or  column  to  column  as  the  case  may  be. 

The  side  walls  are  assumed  to  be  lines  of  constant-z,  the 
back  wall  composed  of  contant-r  and  constant-z  segments  with,  however, 
the  boundary  condition  always  written  only  for  the  constant-r  seg¬ 
ment.  In  this  way,  an  arbitrary  back  wall  can  be  approximated  as 
shown  in  the  following  sketch: 


The  boundary  condition  which  is  applied  at  each  intersection 
of  the  constant-x  lines  with  the  side  walls  and  at  each  intersection 
of  the  constant-z  lines  with  the  back  wall  is 

(Vp 

a  +  bT  =  c  (32) 

where  n  is  the  outward  normal  direction  and  a,  b,  and  c  are 
specified  constants. 

The  energy  input  to  the  ablating  surface  is  determined  from 
a  general  convection  condition.  If  the  mass-transfer  coefficients 
for  all  species  are  equal  and  the  Lewis-Semenov  number  is  unity  we 
have1 


XA  thorough  discussion  of  the  rationale  leading  to  Equation  (2)  is 
given  by  Kendall  and  Rindal,  Reference  4. 
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q  ,  =  p  u  C„  H  +  B 
^net  re  e  H  [_  r 


H  ~  (1  +  B  )H  .  -  oe  T 
a  wj  w 


+  q  ^  +  q_ 

rad  P 


(33) 


where 


r 

H 


€ 

qrad 


heat-transfer  coefficient 
recovery  enthalpy 

enthalpy  of  virgin  ablator  at  wall  temperature 

enthalpy  of  gas  mixture  at  the  wall 

emissivity  of  the  surface 
incident  radiation  heat  flux 

incident  heat  flux  due  to  particle  impingement 
wall  temperature 


The  dimensionless  erosion  rate,  B1,  and  the  wall  gas  enthalpy, 

Hw,  are  obtained  from  an  analysis  of  the  surface  chemical  reactions 
as  described  in  Section  2.2.2.  The  quantities  PeueCH>  Hr»  qrad' 
and  q^  are  properties  of  the  exhaust  stream  and  must  be  derived 
from  an  analysis  of  it.  The  quantities  H  and  e  are  properties 

a 

of  the  wall  material. 

In  order  to  allow  for  some  generality  in  the  problems  which 

can  be  treated,  the  one-moving  boundary  program  is  equipped  with 

three  different  options  for  the  specification  of  surface  boundary 

conditions.  Option  1  relates  specifically  to  the  erosion  of  a 

graphite  wall  in  an  H,  C,  N,  O,  Cl  exhaust.  The  ablating  surface 

heat  flux  is  calculated  by  Equation  (33) .  A  subroutine  is  provided 

which  calculates  B',  H  ,  and  *  :rom  the  set  of  chemical  and 

w  j 


-21- 


traneport  equations  described  in  Section  2.2.2.  The  three  reactions 
given  by  Equations  (2),  (20),  and  (25)  are  presumed  to  proceed  at 
finite  rates. 

The  user  of  the  program  can  either  accept  the  builtrin  reac¬ 
tion  constants  given  in  Section  2.2.2  or  can  specify  his  own  reac¬ 
tion  rates  ir*  the  general  Arrhenius  form 


(34) 


The  specification  of  rate  constants  is  made  by  inputing  the  values 
of  A  and  B  for  each  of  the  three  reactions.  Equations  (2),  (20), 
and  (25) .  Under  this  specification,  reaction  (25)  is  presumed  to 
be  ideal,  that  is,  the  denominator  in  Equation  (29)  is  unity.  The 
builtin  reaction  constants  will  be  revised  throughout  the  remainder 
of  the  current  project  as  (or  if)  newer  or  better  data  become  avail¬ 
able. 

The  remainder  of  the  information  needed  to  evaluate  Equa¬ 
tion  (33)  is  supplied  to  the  program  by  the  way  of  tables  as  func- 
.  tions  of  firing  time.  The  surface  emissivity  and,  of  course,  the 
r-  and  z-direction  thermal  conductivities  are  to  be  provided  as 
*  tabular  functions  of  temperature. 

The  program  modifies  the  heat-transfer  coefficient  for  wall 
blowing  by  the  empirical  relation 


M  molecular  weight  of  the  exhaust  stream  gas  mixture 

M  molecular  weight  of  the  gases  added  at  the  wall 


f 

i 

i 
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The  user  provides  a  table  of  p  u  C„  which  may  be  obtained  by  the 

6  G  xl 

O 

ARM  method  as  discussed  in  Reference  2. 

Under  option  (2),  Equation  (33)  is  still  used  to  calculate 
the  surface  heat  flux,  but  B',  H ,  and  H  can  be  provided  as 
tabular  functions  of  wall  temperature.  In  this  way,  the  same 
moving -boundary  conduction  package  can  be  used  for  materials  or 
exhausts  other  than  those  covered  by  the  chemical  erosion  subroutine. 

Option  (3)  provides  for  specification  of  both  B*  and  as 

tabular  functions  of  firing  time.  This  option  will  be  principally 
used  in  evaluating  test  data  where  the  requisite  information  may  be 
available  from  independent  measurements. 

Complete  details  on  the  use  of  the  program  and  the  various 
boundary-condition  options  will  be  covered  in  a  users  ’  manual.  The 
principal  effort  on  the  one-moving-boundary  program  during  the  next 
quarter  will  be  in  making  minor  modifications  and  the  preparation 
6f  the  users*  manual. 

2.3  Thermal  Response  of  a  Charring  Ablator 
2.3.1  Introduction 

When  a  material,  such  as  a  resin-impregnated  composite,  can 
experience  gross  physical  or  chemical  change  in  the  interior  as  well 
as  at  the  surface,  the  description  of  the  internal  thermal  response 
of  these  materials  to  a  heating  environment  can  no  longer  be  given 
in  terms  of  simple  heat  conduction.  In  addition  to  the  transport 
of  energy  by  conduction,  one  must  consider  the  effects  of  variable 
density,  consumption  of  energy  due  to  formation  of  products  of 
pyrolysis,  and  the  transport  of  energy  through  a  porous  char  material 
by  these  products  of  pyrolysis. 

In  addition  to  the  modification  of  the  equation  describing 
the  internal  thermal  response,  an  analysis  of  these  materials 
requires  a  consideration  of  two  other  phenomena.  First,  the  mecha¬ 
nism  of  the  pyrolysis  itself  must  be  understood  and  quantified. 
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Frequently,  this  requires  a  kinetic  description  of  the  decomposition, 
rather  than  the  assumption  of  a  discontinuous  process  which  occurs 
at  some  threshold  temperature  (ac,  for  instance,  melting) .  Second, 
the  chemical  reactions  between  the  pyrolysis  gases  and  the  gases  at 
the  exposed  surface  must  be  considered. 

All  of  these  phenomena,  collectively,  comprise  the  “response" 
of  a  charring  ablator  to  the  exhaust  environmeht.  All  of  them  and 
the  coupling  among  them  must  be  considered  in  order  to  predict  that 
response.  In  the  following  sections,  the  equations  describing  the 
internal  thermal  response  of  a  charring  ablator  will  be  developed, 
and  the  calculative  technique  used  for  prediction  of  chemical  erosion 
of  the  surface  when  all  the  chemical  reactions  are  fast  (i.e. , 
equilibrium  is  obtained)  will  be  presented.  The  use  of  these  equa¬ 
tions  and  calculative  techniques  to  predict  the  response  of  a  charring 
ablator  to  a  rocket-engine  exhaust  is  the  principal  task  of  the  next 
quarter. 

2.3.2  The  differential  equation  for  the  temperature  field  in  a 
charring  ablator 

Because  the  types  of  materials  comprising  the  class  referred 
to  as  "charring  ablators"  have,  in  general,  a  low  thermal  conduc¬ 
tivity  and  are  used  in  regions  of  a  rocket  nozzle  where  the  axial 
gradients  of  stream  conditions  are  not  severe,  a  one-dimensional 
axisymmetric  body  was  chosen  as  the  representative  shape. 

Consider  the  body  represented  by  the  following  sketch 
(assuming  a  unit  dimension  in  the  z-direction) : 


-  i 


I 

i 
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Let  us  make  an  energy  balance  on  an  annular  ring  of  differential 
thickness,  dr,  as  shown  in  the  sketch.  The  net  energy  transport 
into  the  element  by  conduction  is 


where  K  must  be  throught  of  as  a  local  effective  conductivity. 

If  the  local  rate  of  transpiration  of  pyrolysis  gases  is  denoted 
by  ihg  and  the  enthalpy  of  the  gases  by  Hg,  then  the  net  energy 
transport  into  the  element  by  the  moving  gases  is 

‘  37  <*g  V  dr 
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The  rat®  of  change  of  thermal  energy  storage  in  the  differential 
ring  is 


where  p  and  H  are  the  density  and  enthalpy,  respectively,  of 
the  material  contained  in  the  differential  ring.  Then,  the  equa¬ 
tion  for  the  conservation  of  energy  is  - 

■  He-  -  t  0 « S5)  -  h  »*> 

In  order  to  use  Equation  (34)  for  computation,  we  must  be  able 
to  evaluate  the  material  thermal  and  transport  properties,  pH,  k, 
and  Hg,  as  well  as  the  rate  of  gas  transpiration,  as  functions 
of  position  and  local  temperature.  For  this  purpose,  let  us  con¬ 
sider  each  of  the  terms  in  Equation  (34)  in  greater  detail. 

At  any  instant  in  time,  the  ablating  material  will  be  con¬ 
sidered  as  being  composed  of  two  parts:  the  first,  virgin  plastic 
and  the  second,  completely  decomposed  char.  Notice  that  both  the 
virgin  plastic  and  the  charred  material  may  themselves  be  composites 
of  some  sort.  If  we  denote  properties  of  the  plastic  by  the  sub¬ 
script  p  and  those  of  the  char  by  the  subscript  c,  we  may  write 

P  =  €PPP  +  (1  “  ep)pc  (35) 

where  cp  is  the  volume  fraction  of  undecomposed  plastic  existing 
at  any  instant  in  time.  For  undecomposed  material  ep  is  1,  for 
the  char  layer  ep  is  0,  and  for  other  locations  it  may  be  anywhere 
in  between.  Further,  we  have 


PH  ’  W  Hp  +  (1  -  V  Pc  Hc 


(36) 
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where  and  Hc  are  the  enthalpies  of  the  plastic  and  the  char, 

respectively.  In  order  to  account  for  the  energy  associated  with 
the  process  of  pyrolysis  we  roust  have 


and 


(37) 


(38) 


Implicit  in  the  definition  of  char  is  the  fact  that  both  p  and 

c 

H  must  include  the  gaseous  products  of  pyrolysis  present  in  the 
matrix.  (This  is  a  small  contribution  to  be  sure,  but  the  purist 
will  insist  on  its  inclusion.)  Then,  using  Equation  (36),  we  can 
work  out  the  first  term  in  Equation  (34)  as 


r 


r 


f 


u 


P_€  c 

p  p  p 


+  (1 


€p> 


c  p 

crc 


3 


+ 


(pH 
P  P 


-  PcHc> 
Pc 


(39) 


In  deriving  Equation  (39) ,  p  and  pc  have  been  assumed  to  be 
constants. 

In  order  to  evaluate  the  last  term  in  Equation  (34),  con¬ 
sider  an  enlarged  view  of  the  region  denoted  as  Section  I  in  the 
sketch  on  page  24. 
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We  i«ek  an  expression  for  the  mass  rate  of  flow  of  pyrolysis  gas 
through  any  shell  at  r.  Assume  for  the  moment  that  pyrolysis 
can  be  occurring  anywhere  throughout  the  section  and  that  the 
outer  radius,  r^,  is  at  an  impermeable  surface.  The  mass  of  gas 
crossing  a  cylindrical  surface  a  distance  A^  from  the  back  wall 
is 


r=rb"Ari 


“  2T(rb  -  *1>*1  (It) 


Similarly,  that  crossing  the  surface  a  distance  Arx  +  ir  from 
the  back  wall  is 


r*rb-(Ar1+Ar2) 


r=r, -Ar 
d  x 


+  *•>]  ^  (If) 


If  we  proceed  in  this  way,  adding  up  the  contributions  of  the  annular 
rings  Ar  thick,  until  we  reach  the  surface  at  r,  and  then  allow¬ 
ing  Ar  0  we  have 
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Ag(r) 


(40) 


Note  that  dp/d$  is  itself  a  function  of  r.  We  are  now  prepared 
to  write 


_1_  Hg>  _  a 

2tt  or”  ^ r 


(41) 


The  term  representing  conduction  in  Equation  (34)  is  straight¬ 
forward,  provided  one  can  define  the  conductivi  :y,  k,  for  the  corn- 

► 

posite  material.  If  the  same  gradient  of  temperature  is  assumed  to 
exist,  locally,  in  all  the  component  materials  of  the  composite,  it 
is  easy  to  show  that  the  desired  quantity  is 


k 


I,N 

j 


(42) 


where  €j  is  the  volume  fraction  of  material  j,  and  k^  is  the 
thermal  conductivity  of  material  j.  The  assumption  of  equal  local 
gradients  is  only  strictly  true  for  steady  state, if  the  k^’s  are 
unequal.  Better  results  would  probably  be  obtained  by  using  experi¬ 
mental  data,  when  such  data  are  available. 

If  all  the  terms  appearing  in  Equation  (34)  are  asserobled 
and  the  result  solved  for  dT/d@,  there  results 
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where 


C  =  p  e  C  +  { 1  -  €_)  p  c 

t  P  P  P  P  c 


and 


-  PcHc> 
Pc 


(44) 


(45) 


Whan  ablation  is  occurring  at  the  inner  surface,  the  solution 
of  Equation  (43)  by  finite  difference  techniques  is  facilitated  by 
transforming  to  a  moving  coordinate  system  fastened  to  the  ablating 
surface* 2  Taking 


x  =  r 


G’  =  Q 


(46) 


and  transforming  Equation  (43) ,  there  results 


Notice  that  the  time  derivative  of  temperature  is  taken  at  con¬ 
stant  x,  while  the  time  derivative  of  density  is  at  constant  r. 
The  reason  for  retaining  the  density  term  in  this  form  will  be  dis 
cussed  in  Section  2.3.3 


aThe  roovtives  and  considerations  involved  are  di 3cussed  in  Ref¬ 
erence  1. 
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2.3.3  Thermophysical  properties  of  a  charring  ablator 

In  order  to  effect  a  solution  .to  Equation  (47),  thermo¬ 
physical  data  of  three  different  types  are  required.  These  are: 

(a)  Effective  thermal  conductivities. 

(b)  Density  changes  and  decomposition  rates. 

(c)  Enthalpy-temperature  relationships. 

The  conceptual  difficulty  with  the  apparent  conductivity  of 
a  composite  has  already  been  mentioned.  The  most  satisfactory 
determination  of  conductivity  is  from  experiment,  but  this  is  a 
relatively  difficult  measurement  to  make.  If  the  conductivities 
of  the  materials  comprising  the  composite  are  strongly  tempera¬ 
ture  dependent,  the  effective  conductivity  may  depend  on  the  heating 
rate.  For  most  materials  of  interest,  the  properties  of  the  virgin 
materials  are  fairly  well-established,  and  the  difficulty  is  mainly 
with  partially  or  fully  degraded  material. 

In  general,  the  instantaneous  rate  of  decomposition  of 
phenolic-type  plastics  depends  on  both  the  temperature  and  the 
extent  of  prior  decomposition.  For  example,  Munson  and  Spindler, 
Reference  11,  assume  a  decomposition  law  of  the  form 

=  -  (p  -  pc)n  A  exp(^  If)  (48) 

and  choose  the  constants  to  fit,  experimental  data  on  decomposi¬ 
tion  as  well  as  possible.  In  dealing  with  nylon  cloth-phenolic 
composite,  Kratch,  Kearne,  and  McChesney,  Reference  12,  assume 
that  three  parallel  decompositions  occur,  each  of  which  is  expressed 
in  the  form  of  Equation  (48) .  Using  this,  they  are  able  to  fit  the 
decomposition  rate  data  very  closely. 

It  should  be  noted  that  Equation  (48)  (or  its  counterparts 
in  a  more  elaborate  decomposition  model)  applies  to  a  point  fixed 
in  the  decomposing  material;  that  is,  at  some  point  in  the  ablating 
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body  there  exist  some  instantaneous  values  of  p  and  T  which 
will  uniquely  determine  dp/d&.  Therefore,  the  value  resulting 
from  Equation  (48)  represents  the  (dp/d0)r  term  appearing  in 
Equation  (47) . 

The  enthalpies  and  specific  heats  of  the  virgin  plastic 
and  the  char,  as  functions  of  temperature,  can  either  be  determined 
from  direct  calorimetric  measurements  or  computed  from  data  on  the 
component  molecules  if  such  data  are  available.  The  enthalpy  of 
the  mixture  of  pyrolysis  gases  can  be  calculated,  provided  some 
additional  assumptions  regarding  chemical  and  thermal  state  are 
made.  It  appears  reasonable  to  assume  that  the  gases  achieve 
both  thermal  and  chemical  equilibrium  with  the  char  layer. 

If  the  density  of  the  material  at  any  point  and  at  any 
instant  in  time  is  known,  the  apparent  fraction  of  undecomposed 
plastic,  6  ,  can  be  calculated.  It  should  be  appreciated  that 
€p  is  a  convenient  fiction  that  truly  represents  the  state  of 
affairs  only  at  the  values  0  and  1. .  It  may  be  possible  to  estab¬ 
lish  the  adequacy  of  the  approximation  by  way  of  Equation  (36), 
if  enthalpy  data  are  available  for  partially  degraded  material, 
as  well  as  for  virgin  material  and  char. 

One  area  of  continuing  effort  on  this  program  is  a  search 
for  the  requisite  thermophysical  data  for  the  various  materials 
of  interest.  It  is  hoped  that  analyses  such  as  that  outlined 
previously  will  stimulate  the  acquisition  of  the  necessary  ma¬ 
terial  properties. 

2.  3.4  Surface  chemical  reactions 

In  order  to  establish  a  surface  boundary  condition  for  the 
internal  conduction  program,  it  is  necessary  to  consider  the  requi¬ 
site  inputs  and  resultant  outputs  involved  in  a  calculation  of  the 
surface  state.  At  each  time  step,  certain  information  is  avail¬ 
able  from  the  internal  conduction  solution  which  may  be  used  to 
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evaluate  this  surface  state.  Included  are: 

(a)  The  rate  at  which  pyrolysis  gas  is  being  supplied  to 

the  surface  as  a  consequence  of  in-depth  decomposition 
of  resin  (for  a  charring  ablator) . 

(b)  A  relation  between  the  surface  temperature  and  the 

surface-temperature  derivative;  for  example,  using  the 
temperature  of  a  node  ix  below  the  surface,  this 
relation  is  of  the  general  form 


f(T  ) 
w 


(49) 


In  general,  several  other  required  parameters  are  established 
either  absolutely  or  through  some  Known  or  calculated  functional 
dependence.  These  parameters  include: 

(a)  The  elemental  composition  of  the  char  layer. 

(b)  The  elemental  composition  of  the  pyrolysis  gas. 

(c)  The  elemental  composition  of  the  boundary- layer-edge  gas. 

(d)  ,The  enthalpy  of  the  boundary-layer-edge  gas. 

(e)  The  incident  radiation  flux  upon  the  wall. 

(f)  Heat-  and  mass-transfer  coefficients. 

(g)  The  blowing  correction  factor  on  the  preceding. 

(h)  The  enthalpy  of  the  char  as  a  function  of  its  tempera¬ 

ture. 

(i)  The  enthalpy  of  the  adjacent  pyrolysis  gas  (immediately 

below  the  surface)  as  a  function  of  its  temperature. 

(j)  The  thermal  conductivity  of  the  char  as  a  function  of 

temperature. 

(k)  The  emissivity  of  the  char  as  a  function  of  temperature. 
Some  of  the  above  parameters  may  be  given  explicitly  as  functions  of 
time  without  further  complicating  the  problem. 


\ 
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From  this  information  it  is  necessary  to  establish  the  surface 
.temperature  and  the  char-removal  rate.  The  present  procedure  at 
jvioya  utilizes  mass  and  energy  balances  plus  an  assumption  of  sur¬ 
face  chemical  equilibrium.  In  the  following  paragraphs  the  func¬ 
tional  relations  will  be  presented  briefly,  together  with  some 
eample  results. 

By  performing  an  elemental  mass  balance  on  each  of  the  I 
elements  at  the  surface,  it  is  possible  to  formulate  an  expression 
for  the  mass  fractions  of  each  of  the  elements  contained  in  the  gas 
immediately  adjacent  to  the  material  surface.  These  relations  are 
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where 


+ 

i  -  q  c. 
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PeU  C 
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(51) 


for  the  case  of  no  condensed-phase  deposition  or  removal  from  the 
|ur£ace  and  equal  effective  binary  diffusion  coefficients  within 
the  boundary  layer.  Based  on  the  presumed  input  information,  the 
t  equations  of  this  type  contain  1+1  unknowns,  Kiw  and  ftc« 
The  energy  balance  given  by  Equation  (33)  can  be  written 
somewhat  more  generally  (in  the  absence  of  particles)  as 
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where  the  first  equality  is  a  consequence  of  the  internal  conduction 
solution.  This  equation  introduces  two  new  unknowns,  namely,  the 
enthalpy  of  the  wall  gas,  Hw,  and  the  wall  or  surface  temperature, 

v 

The  assumption  of  surface  equilibrium,  both  among  the  gas 
species  adjacent  to  the  surface  and  between  them  and  the  surface 
material,  introduces  the  final  required  relations.  It  is  well- 
established  that,  given  the  elemental  composition  of  a  gas  and  its 
temperature  and  pressure,  a  definitive  set  of  relations  can  be  formu¬ 
lated  to  establish  the  state  of  the  system.  In  the  present  case 
this  would  imply  that 

Surface  gaseous  state  =  f(K.  ,  T  ,  p)  (53) 

1W  w 


using  the  requirement  that  gas  phase  equilibrium  exists  among  the 
species  at  the  surface.  The  requirement  for  equilibrium  between 
the  gas  phase  and  the  surface  material  has  yet  to  be  included.  The 
introduction  of  this  constraint  demands  the  removal  of  ^ne  of  the 
existing  constraints  in  Equation  (53) .  Consequently,  we  may  write 
for  the  ablation  problem 

Surface  gaseous  state  =  f(K.  ,  p)  (54) 

iw 


or  more  specifically 


=  f(K.  ,  p) 
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(55) 

=  f  (K .  ,  p) 
iw  r 

(56) 

Thus,  a  definitive  set  of  relations  exists  in  the  1+3  unknowns, 

'V 

K.  ,  T  ,  H  ,  and  ft  .  and  the  1+3  Equations  (50),  (52),  (55),  and 
iw  w  w  c 

(56). 


The  solution  of  the  system  of  relations  represented  by 
Equation  (54)  has  been  programed  at  Vidya  for  solution  on  the 
IBM  7C94  computer.  The  program  has  the  capability  of  selecting 
the  appropriate  surface  material  from  a  list  of  possible  candidates. 
For  example,  the  exposed  surface  of  a  silica-phenolic  ablator  can 
be,  under  various  conditions,  either  silicon  carbide,  carbon,  liquid 
silicon,  or  silica.  A  sample  of  the  output  from  this  program  for 
silica  phenolic  is  shown  in  Table  I.  This  sample  computer  output 
was  selected  to  show  the  transition  from  one  surface  material  to 
another  as  the  char-removal  rate  varies  at  fixed  pyrolysis  gas 
efflux.  The  variation  of  surface  temperature  indicates  why  it  has 
been  discarded  as  an  independent  variable  in  these  calculations. 

The  establishment  of  the  surface  boundary  state  has  thus 
been  reduced  to  an  iteration  based  upon  the  selection  of  an  initial 
char-recession  rate,  Ac,  evaluation  of  the  elemental  mass  fraction 
of  the  surface  gases,  K.  ,  with  Equation  (50) ,  the  establishment  of 

1W 

the  wall  temperature  and  enthalpy  through  the  use  of  the  described 
computer  program,  and,  finally,  the  checking  of  these  values  by 
use  of  the  energy  equation  (Eq.  (52)).  The  actual  procedure  for 
this  iteration  is  based  upon  a  simple  Newton  method  approach. 

The  foregoing  discussion  has  shown  how  a  multidependent 
system  has  been  reduced  to  a  simple  iteration  based  on  a  single 
parameter  by  use  of  a  computer  program,  written  at  Vidya,  for  the 
evaluation  of  a  surface  equilibrium  state  of  an  ablating  material. 
Because  of .the  complete  generality  of  the  program,  it  is  applicable 
to  all  material  types,  requiring  only  the  knowledge  of  conventional 
thermodynamic  properties  as  contained,  for  example,  in  Reference  13. 

3.  EXPERIMENTAL  STUDIES 
3. 1  Introduction 

The  choice  of  a  rocket-nozzle  liner  material  requires  a  know¬ 
ledge  of  the  behavior  of  appropriate  materials  in  the  combustion 
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product  environment.  This  knowledge  must  be. supplied  by  both  a 
theoretical  treatment  of  the  materials  response  and  experimental 
materials  testing.  As  indicated  in  Figure  1,  the  theoretical 
treatment  requires  appropriate  input  of  basic  data  for  a  success¬ 
ful  and  accurate  solution.  Some  of  this  input  must  be  determined 
experimentally.  Hence,  in  the  experimental  phases  of  this  program, 
several  areas  in  which  data  are  currently  lacking  are  being  investi¬ 
gated.  These  areas  include  chemical  reactions  and  particle  impact 
as  they  affect  surface  erosion.  The  experimental  effort  is  con¬ 
veniently  divided  into  two  separate  programs: 

(1)  The  chemical-erosion  program;  the  investigation  of 
chemical  effects  in  the  erosibn  process  and  the  development  of  a 
simple  materials-testing  technique. 

(2)  The  particle-impact  program;  the  investigation  of  the 
effects  of  particle  impingement  on  erosion. 

As  indicated  above,  a  secondary  purpose  of  the  chemical- 
erosion  program  is  to  investigate  a  relatively  simple  materials- 
testing  technique.  The  usual  materials-testing  techniques  have 
taken  three  general  forms:  full-scale  rocket  tests,  small-scale 
rocket  tests,  and  splash  tests  on  materials  samples.  The  last 
technique  must  be  regarded  as  only  qualitative  since  only  a  rough 
comparison  between  the  performance  of  different  materials  is 
possible.  With  the  advent  of  large  rocket  motors,  full-scale 
rocket  tests  for  materials  evaluation  become  prohibitive.  Small- 
scale  rocket  tests  do  not  permit  the  flexibility  of  varying  test 
parameters  independently  and  are  also  relatively  expensive.  A 
flexible,  inexpensive  yet  quantitative  experimental  materials- 
evaluation  technique  is  therefore  desirable,  and  its  development 
is  the  subject  of  one  part  of  the  experimental  effort. 

The  presentation  which  follows  covers  the  areas  of  experi¬ 
mental  activity  during  the  third  quarter  of  the  project.  In  the 
chemical-erosion  test  program,  43  test  firings  which  were  a  direct 
part  of  the  program  were  performed  during  the  report  period.  This 
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does  not  include  a  large  number  of  checkout  tests  which  were  also 
run.  These  43  tests  included  the  firing  of  24  nozzles  of  graphite, 
graphite  phenolic,  and  silica  phenolic  and  16  heat-transfer  and 
pressure-distribution  calibration  tests.  This  testing  completed 
the  program  on  Mixtures  4  and  5  (see  Ref.  2) .  Testing  will  be 
continued  for  Mixtures  1-3  later  in  the  next  quarter.  The  program 
will  be  stopped  in  the  early  part  of  the  quarter,  however,  to 
allow  a  preliminary  test  series  in  the  particle-impact  program  as 
discussed  below. 

The  effort  in  the  particle-impact  test  program  consisted 
of  the  definition  and  investigation  of  instrumentation,  particularly 
for  measurement  of  particle  velocity,  and  the  start  of  the  setup 
for  the  program.  The  program  effort  for  the  immediate  fu  are  will 
be  the  setup  for  and  performance  of  a  preliminary  test  series  to 
get  the  "feel"  of  the  program  and  to  define  problem  areas,  if  any. 

If  problems  do  come  up,  the  appropriate  fixes  will  be  made  while 
the  chemical-erosion  program  is  continued. 

3.2  Chemical-Erosion  Test  Program 
3.2.1  Introduction 

The  primary  objective  of  the  chemical-erosion  test  program 
is  the  determination  of  the  effects  of  specific  chemical  reactions 
on  the  erosion  of  rocket  liner  materials.  The  results  obtained, 
in  addition  to  being  of  general  importance,  provide  necessary 
guidance  to  the  theoretical  phase  of  the  subject  contract.  A 
secondary  objective  is  the  definition  of  an  inexpensive  technique 
for  the  quantitative,  as  opposed  to  qualitative,  evaluation  of 
candidate  materials  for  rocket  nozzles.  The  technique  is  expected 
to  be  a  powerful  tool  for  the  material  fabricator  and  the  nozzle 
designer. 
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The  program,  in  brief,  consists  of  an  appropriate  series  of 
tests  on  ablative  materials  utilizing  the  five  gas  mixtures  discussed 
in  References  1  and  2.  These  gas  mixtures  closely  approximate  the 
temperature-enthalpy  variation  of  a  typical  solid -propellant  exhaust 
gas  and  range,  chemically,  from  a  duplication  of  the  actual  combus¬ 
tion  products  to  a  mixture  of  inert  gases.  The  Vidya  arc-plasma 
generator  supplies  the  energy  input  to  the  gas  mixtures  for  simu¬ 
lation  of  the  rocket  environment.  The  test  materials  are  graphite, 
silica  phenolic,  and  graphite  phenolic.  An  axisymmetric-nozzle 
configuration  is  used  for  all  tests.  The  heat  transfer  and  pressure 
distributions  for  a  nonablating  wall  of  the  same  geometry  as  the 
test  nozzle  are  determined  with  sets  of  heat-transfer  calibration 
nozzles  and  pressure-distribution  calibration  nozzles. 

Testing  during  the  report  period  included  24  nozzle  firings 
of  graphite,  graphite  phenolic,  and  silica  phenolic.  Several  heat- 
transfer  and  pressure-distribution  calibration  tests  were  also  per¬ 
formed.  The  results  of  these  tests  are  presented  in  the  following 
sections.  In  addition  to  this  effort  which  was  a  direct  part  of 
the  chemical-erosion  program,  several  checkout  teats  including 
nine  graphite-nozzle  firings  were  performed  to  eliminate  a  swirl 
component  in  the  exit  flow  and  to  establish  a  high-confidence  level 
in  the  test  results.  Because  of  unanticipated  problems  in  elimi¬ 
nating  the  swirl,  the  anticipated  completion  of  the  chemical-erosion 
program  by  the  end  of  this  report  period  was  not  achieved.  Approxi¬ 
mately  half  of  the  program  was  completed,  however;  all  tests  with 
Mixtures  4  and  5  were  performed.  The  chemical-erosion  program 
effort  will  be  reduced  in  the  immediate  future  to  allow  an  initial 
testing  effort  in  the  particle-impact  test  program.  This  change 
in  emphasis  is  being  made  in  order  to  get  a  feel  for  the  problems, 
conduct,  and  nature  of  the  results  in  the  particle-impact  program. 
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3.2.2  Facility  modifications 

Two  significant  facility  modifications  were  made  during  the 
report  period,  one  to  eliminate  the  swirl  component  in  the  exit  flow 
of  the  plasma  generator  and  the  other  to  improve  the  power  supply 
capability  of  the  facility.  These  modifications  and  additions  are 
discussed  below. 

As  presented  in  Reference  1,  the  discussion  of  the  arc-plasma 
generator,  arc  stability  and  the  integrity  of  the  electrodes  is 
maintained  by  inducing  arc  rotation  with  a  solenoidal  magnetic  field. 
This  arc  rotation  induces  an  undesirable  angular  velocity  or  swirl 
component  in  the  flow  that,  under  some  conditions,  can  be  quite 
significant.  Tests  performed  early  in  the  experimental  program 
indicated  that  the  swirl  could  be  kept  to  an  apparently  acceptable 
level  by  introducing  the  primary  gas  tangentially  in  a  direction 
counter  to  that  induced  by  arc  rotation.  These  teste  were  performed 
on  flat-faced  nylon  disks  placed  normal  to  the  plasma  generator  exit 
flow.  The  flow  pattern  indicated  by  the  melted-nylon  surface  showed 
that  the  swirl  could  be  brought  down  to  a  negligible  level  by  using 
sonic  or  close  to  sonic  jets  for  introducing  the  primary  gas.  This 
splash  flow  apparently  was  not  a  critical  indicator  of  the  presence 
of  swirl,  however,  since  close  inspection  of  the  throat  entry  of 
the  first  nozzle  fired  (Test  No.  795,  reported  in  R$f.  2)  revealed 
a  pattern  that  indicated  the  presence  of  a  significant  swirl  com¬ 
ponent. 

Further  testing  on  eliminating  the  swirl  was  therefore  begun 
early  in  this  report  period.  The  magnetic-field  strength  (which 
dictated  the  magnitude  of  the  swirl  component)  was  first  reduced. 

At  a  field  strength  one-fourth  of  that  usually  used,  the  swirl  as 
indicated  in  nozzle  firings  finally  became  negligible.  Due  to  the 
low  field  strength,  however,  the  plasma-generator  operation  was 
somewhat  unstable  -  sufficienctly  so  that  operation  at  these  con¬ 
ditions  was  unacceptable.  Tangential  injection  of  the  secondary 


gas-flow  was  then  attempted  with  successful  results.  By  injecting 
the  secondary  gas  in  a  direction  counter  to  the  arc-induced  rota¬ 
tion  of  the  primary  gas,  the  swirl  was  eliminated.  No  stability 
problems  existed,  since  swirl  elimination  was  done  at  the  standard 
magnetic-field  strength  of  6000  gauss.  The  facility  modification 
required  to  accomplish  this  was  simply  the  fabrication  of  a  new 
injection  ring.  The  characteristics  of  the  new  ring  are  as  follows: 
the  angle  of  injection,  the  angle  between  a  radial  line  from  the 
plasma-generator  centerline  to  the  injection  point  and  the  direction 
of  secondary  gas-injection,  is  42°;  16  injection  ports  are  used; 
and  the  injection  jets  are  changeable  to  allow  swirl  elimination 
at  different  secondary-gas  conditions. 

The  increase  in  power  supply  capability  mentioned  above  was 
the  addition  of  a  diesel  generator  to  the  power  supply  system.  This 
was  an  in-house  sponsored  effort  as  part  of  Vidya’s  continuing 
program  to  increase  the  capabilities  of  the  arc-plasma  generator 
facility.  The  diesel  unit  was  put  into  operation  late  in  the 
report  period  and  is  capable  of  producing  up  to  500  kw  of  electri¬ 
cal  power.  A  tandem  hookup  with  the  battery  power-supply  is  used 
to  eliminate  the  power  droop  inherent  in  the  battery  system  alone. 

By  "topping"  the  batteries  with  the  diesel  (the  batteries  and 
diesel  in  parallel)  constant-power  operation  in  excess  of  one  mega¬ 
watt  is  possible  for  extended  operating  times.  The  frequency  at 
which  tests  can  be  performed  has  also  been  increased  since  the 
diesel  unit  is  also  used  to  recharge  the  battery  supply  at  a  much 
higher  rate  than  previously  possible.  This  new  power-supply  capa¬ 
bility  has  already  been  used  in  the  chemical-erosion  test  program 
for  the  intermediate-enthalpy  test  firings,  the  last  firings  that 
were  performed  during  the  quarterly  report  period. 

3.2.3  Instrumentation  and  data  reduction 

The  instrumentation  and  data-reduction  techniques  have  been 
presented  for  the  most  part  in  References  1  and  2;  these  will  not 
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be  repeated  here.  Three  particular  areas,  the  X-ray  technique  for 
locating  thermocouple  positions  in  the  test  nozzles,  the  data- 
reduction  computer  programs,  and  the  pyrometer  measurement  of  sur¬ 
face  temperature  have  not  been  discussed  previously,  however,  and 
are  presented  below. 

As  discussed  in  Reference  2,  the  thermocouple  instrumentation 
for  measurement  of  the  internal-temperature  distributions  in  the 
test  nozzles  consists  of  very  small  diameter  probes  (0.035  inch). 
These  probes  are  inserted  into  the  test  nozzles  so  that  they  follow 
approximately  an  isotherm  in  the  region  of  the  thermocouple  junction 
in  order  to  minimize  the  isotherm-distortion  effect  on  the  tempera¬ 
ture  readings.  Because  of  the  small  hole  size  and  large  depth, 
drift  of  the  drill  during  the  hole-drilling  operation  can  occur, 
particularly  in  the  phenolics  where  there  are  material  nonuni forroi- 
ties.  The  location  of  the  bottom  of  the  hole  (the  location  of  the 
thermocouple  bead)  may  therefore  not  be  where  expected  from  measure¬ 
ment  of  the  point  of  drill  entry.  X-ray  photographs  of  each  instru¬ 
mented  test  nozzle  are  therefore  taken  at  a  commercial  inspection 
laboratory.  The  technique  for  determining  the  thermocouple  loca¬ 
tions  from  the  X-ray  photographs  is  illustrated  in  the  following 
sketch 4 


Film 

plane 
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Correction  roust  be  roade  for  the  effect  of  the  film  plane  and  the 
thermocouple  plane  being  non-coplanar . -  the  thermocouples  being 
above  the  film  plane.  In  order  to  make  this  correction,  however, 
the  following  requirements  must  be  met:  the  X-ray  source  must 
approximate  a  point,  the  film  plane  and  the  nozzle  must  be  level, 
and  the  nozzle  centerline  must  coincide  with  the  normal  from  the 
"film  plane  to  the  source.  Unfortunately,  commerical  (as  well  as 
medical)  X-ray  equipment  is  not  ordinarily  set  up  within  these 
requirements.  A  number  of  tests,  measurements,  and  adjustments 
were  therefore  made  on  the  equipment  to  allow  the  required  pre¬ 
cision  measurements.  The  X-ray  unit  used  has  a  small-source 
diameter,  0.5  mm.  The  distance  from  the  source  to  the  film  plane 
was  measured,  and  is  48  inches.  The  film  plane  was  leveled  and 
the  "bulls-eye"  on  the  film  plane  was  precisely  located.  The 
required  capability  for  accurately  determining  the  thermocouple 
locations  was  therefore  established. 

A  typical  X-ray  print  is  shown  in  Figure  4.  The  three  dark 
lines  which  terminate  approximately  perpendicular  to  a  radial  line 
from  the  nozzle  center  are  wires  inserted  in  the  probe  holes  to 
simulate  the  thermocouple  probes.  The  dark  rectangular  patches 
are  calibration  shims  of  known  length  for  length  calibration  and 
for  a  check  on  the  accuracy  of  the  calculated  locations.  One  of 
the  long  shims  is  on  the  top  surface  of  the  nozzle  and  the  other 
is  on  the  film  plane.  By  making  the  appropriate  correction  on 
the  measured  length  of  the  top  shim,  the  accuracy  of  the  method 
is  checked  for  each  nozzle.  A  cylindrical  plug  is  inserted  into 
the  tubular  throat  region  to  define  accurately  its  location;  with¬ 
out  it,  the  definition  of  the  throat  surface  is  somewhat  fuzzy. 

For  actual  measurement,  the  X-ray  negative  is  enlarged 
approximately  three  times.  The  film-plane  shim  is  then  used  for 
a  length  calibration.  The  error  in  the  calculated  distance  from 
the  nozzle  centerline  to  the  probe  tip  is  felt  to  be  within 
+3  percent. 
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In  order  to  speed  up  and  simplify  the  data-reduction  process 
for  the  nozzle  firings  and  calibration  tests >  a  number  of  computer- 
program  additions  and  modifications  were  written  during  the  report 
period  for  use  on  the  Vidya  IBM  1620  computer.  A  heat-flux  calcu¬ 
lation  routine  was  added  to  the  standard  arc-plasma-generator  data- 
reduction  program  for  use  in  the  reduction  of  the  heat-transfer 
calibration  test  data.  The  average  heat  flux  for  each  calibration 
nozzle  segment  is  calculated  and  read  out  along  with  the  standard 
arc-plasma-generator  operating  data. 

A  new  program  was  written  to  check  nozzle  performance,  to 

calculate  instantaneous  surface  recession,  and  to  calculate  the 

instantaneous  throat  heat-  and  mass-transfer  coefficients.  The 

nozzle  performance  is  checked  by  a  calculation  of  the  nozzle 

coefficient,  C  ,  (see  Eq.  (63) ,  Ref.  2)  for  each  test.  Values  not 
n 

in  an  appropriate  range  indicate  some  problem  in  the  test  or  in 
the  data  reduction  (the  values  of  Cn  typically  fall  in  the  range 
0.90  to  0.98).  The  instantaneous  surface  recession  is  calculated 
from  the  plenum  pressure  at  a  number  of  discrete  times  during  a 
test  by  first  calculating  the  throat  area  at  each  time,  then  the 
throat  diameter  obtained  therefrom  and,  from  this  and  the  initial 
throat  diameter,  the  surface  recession  (see  pp  57  and  58,  Ref.  2) . 
The  heat-transfer  coefficient  is  calculated  from  the  simplified 
Bartz  equation  (see  Section  3. 2.4.3)  at  each  time  for  which  sur¬ 
face  recession  is  determined.  The  mass-transfer  coefficient  ij 
calculated  from  the  analogy  between  heat  and  mass  transfer  (see 
Section  3. 2. 4. 4. 4).  The  required  transport  properties  are  pro¬ 
gram  inputs  along  with  the  thermodynamic  data  obtained  from  the 
arc-plasma-generator  operation  conditions.  This  program  is  also 
used  as  a  check  on  the  nozzle  coefficient  in  the  pressure  distri¬ 
bution  and  heat-transfer  calibration  tests,  and  for  calculation 
of  heat-transfer  coefficient  in  the  heat-transfer  tests. 

An  existing  curve-fit  program  was  modified  for  calculation 
of  the  instantaneous  surface-recession  rate.  The  calculated 
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values  of  surface  recession  are  input  to  the  program  and  a  least 
squares  fit  of  these  points  made.  The  slope  of  the  fitted  curve, 
the  surface-recession  rate,  is  calculated  and  read  out  together 
with  the  curve  fit  itself. 

A  computer  program  was  also  written  to  reduce  the  pre-  and 
post-test  nozzle  profile  data.  The  axial  variation  of  nozzle 
diameter  before  and  after  test  and  the  axial  and  circumferential 
variation  of  surface  recession  are  calculated.  The  program  input 
is  the  two  sets  of  readings  (before  and  after  test)  from  four  axial 
traverses,  each  90°  apart  circumferentailly,  of  a  pivot  arm-dial 
gauge  micrometer. 

A  two-color  optical  pyrometer  is  being  used  in  all  nozzle 
firings  in  an  attempt  to  get  a  direct  indication  of  wall  tempera¬ 
ture.  The  pyrometer  is  set  so  that  it  looks  into  the  nozzle  throat 
at  the  throat  wall.  The  two-color  principle  has  the  advantage  of 
requiring  no  emissivity  correction.  The  pyrometer  senses  the 
intensity  at  two  discrete  wavelengths  and,  from  the  ratio  of  the 
intensities,  identifies  the  temperature  corresponding  to  the 
black-body  curve  that  has  this  intensity  ratio.  Since  it  senses 
temperature  from  a  ratio  of  intensities  rather  than  from  an  abso¬ 
lute  intensity,  the  temperature  readout  is  independent  of  emis¬ 
sivity  if  the  object  is  radiating  as  a  gray  body.  This  type  of 
pyrometer  also  has  its  disadvantages,  however,  the  most  serious 
being  its  sensitivity  to  stray  radiation.  Any  stray  radiation 
intercepted  by  the  pyrometer  will  generally  have  different  rela¬ 
tive  intensities  than  the  object  of  interest  and  therefore  will 
throw  the  pyrometer  indication  off.  Because  of  this  effect,  the 
validity  of  the  pyrometer  data  obtained  by  looking  through  the 
luminous  arc-plasma  stream  is  somewhat  questionable. 3 


3Thermocouple  instrumentation  is,  of  course,  also  being  used  for 
determination  of  wall  temperature.  The  pyrometer  complements 
these  measurements. 
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3.2.4  Test  results  and  discussion 
3.2.4. 1  Introduction 

The  test  program  on  Mixtures  4  and  5  was  completed  during 
the  report  period  and  the  majority  of  results  are  presented  herein. 
Because  of  incomplete  data  reduction  on  some  tests  at  the  time  of 
writing,  some  results  are  necessarily  omitted.  Twenty-four  test 
nozzles  of  graphite,  graphite  phenolic,  and  silica  phenolic  were 
fired  as  a  direct  part  of  the  program.  Appropriate  heat-transfer 
and  pressure-distribution  calibration  tests  were  also  performed. 

The  complete  test  program  is  presented  in  Table  II.  This  table 
does  not  include  a  number  of  other  tests,  including  nine  graphite- 
nozzle  firings,  which  were  performed  during  the  report  period  in 
preparation  and  checkout  for  the  actual  program.  Note  that  the 
test  program  was  altered  somewhat  from  that  originally  presented 
in  Reference  2.  This  was  because  Graphitite  GX  exhibited  no  sur¬ 
face  recession  at  all  with  Mixture  5  (all  inert  gases)  at  the  high- 
temperature  condition,  and  therefore  the  two  graphite-nozzle 
firings  at  less  severe  conditions  were  eliminated;  two  firings 
with  ATJ  graphite.  Test  Nos.  1000  and  1005,  were  substituted  in 
their  place  to  obtain  a  comparison  between  ATJ  and  Graphitite  GX. 

The  test  results  are  presented  and  discussed  in  the  follow 
ing  sections.  The  experimental  pressure  distributions  are  pre¬ 
sented  in  terms  of  the  pressure  ratio,  p/pQ,  and  are  compared  with 
one-dimensional  isentropic  theory.  The  experimental  heat-transfer 
results  are  presented  in  the  form  of  heat-transfer  coefficients 
and  are  compared  with  the  axial  distribution  predicted  by  the 
simplified  Bartz  method.  The  results  of  the  nozzle  firings  are 
presented  and  discussed  in  detail.  These  results  include  post¬ 
test  observations  of  surface  condition,  char  depths,  and  surface 
recession,  the  time  variation  of  surface  recession  and  surface- 
recession  rate,  and  internal  and  wall-surface  temperature  uistoriee 
Wherever  possible,  comparisons  are  made  with  theoretical  predic¬ 
tions. 
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3. 2. 4. 2  Pressure-distribution  tests 

Pressure-distribution  calibration  tests  were  run  for  the 

D  0 

complete  range  of  test  conditions,  nominally  2000  K  to  3500  K, 
with  both  Mixtures  4  and  5.  The  description  of  the  nozzles  and 
conduct  of  tests  is  presented  in  References  1  and  2  and  is  not 
repeated  here.  The  test  results  are  presented  in  figures  5  and  6 
and  are  discussed  below.  The  test  conditions  are  indicated  in 
each  figure-  as  is  the  theoretical  one-dimensional  isentropic 
pressure  distribution  for  each  nozzle  contour. 

As  seen  from  the  figures,  the  experimental  results  arc 
characterized  by  satisfactory  agreement  with  one-dimensional 
isentropic  theory  in  the  convergent  entry  section  and  the  divergent 
exit  section  of  the  nozzles  (except  as  noted  below)  but,  as 
expected,  substantial  disagreement,  particularly  for  the  0. 3-inch- 
diameter  throat  nozzle,  is  apparent  in  the  tubular  throat.  As 
'Hscussed  in  detail  in  Reference  2,  this  discrepancy  is  due  to 
machining  nonuniformities  and  boundary-layer  growth.  The  pre¬ 
dicted  pressure  drops  based  on  c  nsideratioi  of  these  effects 
were  presented  in  Reference  2  for  the  0.3-and  0.4-inch-diameter 
throat  nozzles  and  was  calculated  during  this  report  period  for 
the  0.5-inch  throat  nozzle.  A  comparison  of  the  predicted  values 
and  the  experimental  results  is  presented  in  the  table  below: 


D* 

(in.) 

Predicted  pressure  drop 
(psi) 

Measured  pressure  drop 
(psi) 

0.3 

52 

49 (956) ' 

55 (955) 

.4 

21 

10(954}' 

11  (953)  ' 

10 (1030) ' 

11 (1029) 

.5 

9 

6  (952)  ’ 

6  (951) 

Numbers  in  parentheses  indicate  test  numbers. 
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Discrepancies  between  measured  and  calculated  values  are  felt  to 
be  due  to  the  simplified  approach  used  for  predicting  the  pressure 
drop  (Ref.  2) ,  and  to  the  small  but  progressive  copper  oxide  build¬ 
up  which  voccurs  during  each  test  causing  a  slight  change  in  the 
nozzle-throat  profile  and  therefore  in  the  pressure  drop. 

Relatively  poor  agreement  between  experiment  and  theory 
is  evident  in  the  divergent  exit  section  of  the  high-temperature 
condition  tests  on  the  0.3-  and  0.4-inch-diameter  throat  nozzles, 
Figures  5U'  and  {b)  and  6(a)  and  (b)  .  No  explanation  for  this 
discrepancy  is  apparent.  The  possibility  of  z  calibration  shift 
in  the  pressure  transducer,  a  leak  in  the  pressure  line,  and  an 
obstruction  in  the  line  were  all  investigated  with  a  negative 
result.  Also,  the  0.5-inch-throat  nozzle  run  just  prior  to  these 
tests  (Figs.  5(c)  and  6(c))  demonstrates  favorable  agreement  as 
does  the  0.4-inch  nozzle  run  at  the  low-temperature  condition  later 
in  the  program  (Figs.  5(d)  and  6(d)). 

Results  at  four  of  the  nominal  test  conditions  presented 
herein.  Figures  5(a)  and  (b)  and  6(a)  and  (b) ,  were  also  presented 
in  Reference  2.  These  four  tests  were  repeated,  however,  since 
those  tests  presented  in  Ref  '.  *ence  2  were  run  before  the  swirl 
component  in  the  flow  was  eliminated.  From  s  comparison  of  results, 
however,  the  swirl  component  was  found  to  have  had  little  effect 
on  the  pressure  distribution. 

3. 2.4. 3  Heat-transfer  tes:s 

As  discussed  in  detail  in  Reference  2,  the  heat-transfer 
calibration  nozzles  consist  of  stacked  calorimeter  segments  which 
form  the  desired  nozzle  contour.  The  average  steady-state  heat 
flux  for  each  segment  is  determined  from  the  known  heated  wail 
surface  area  of  the  segment  and  the  heat  rate  into  the  segment 
calculated  from  the  measured  cooling-water  flow  rate  and  tempera¬ 
ture  rise.  The  reader  is  referred  to  Reference  2  for  further 
details  of  the  nozzle  construction  and  data  acquisition. 
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Tests  weio  run  on  the  0.3-,  0.4-,  and  0. 5-ineh-diameter 
throat  nozzles  with  Mixtures  4  and  5  over  the  complete  free-stream 
temperature  range  of  the  program,  nominally,  2000°  K  to  3500J  K, 

The  experimental  results  for  the  high-temperature  condition  are 
presented  in  Figure  7  and  8.  The  data  for  the  low -temperature 
tests,  which  were  obtained  on  the  0.4-inch-diameter  throat  nozrie 
only,  are  in  the  data-reducfcion  process  and  are  not  yet  available 
for  presentation  herexn.  In  each  figure,  the  nozzle  contour  and 
segment  positions  are  indicated  for  reference.  The  experimental 
results  are  presented  in  terms  of  the  heat-transfer  coefficient, 
defined  as  q/(Hr  -  H^) ,  where  the  recovery  enthalpy  is  approxi¬ 
mated  by  the  stagnation  enthalpy  (a  reasonable  approximation  here). 
The  wall  enthalpy  was  calculated  from  the  measured  heat  flux  by 
considering  the  water  side  heat  transfer  and  the  heat  transfer 
through  the  copper  wall.  In  the  variable  area  segments,  the 
convergent  and  divergent  sections,  the  experimental  heat-transfer 
coefficients  are  plotted  at  the  axial  point  where  the  predicted 
average  heat-transfer  coefficient  for  the  segment  is  equal  to  the 
predicted  local  coefficient  based  on  the  simplified  Bartz  tech¬ 
nique  discussed  below. 

The  predicted  axial  variation  of  heat -transfer  coefficient 
is  also  included  in  each  plot  for  reference.  This  variation  was 
calculated  from  the  simplified  Bartz  expression  (Refs.  14  and  2) , 
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where  the  prime  superscripts  indicate  evaluation  of  properties  at 
the  film  enthalpy4. 


The  reference  enthalpy,  as  presented  in  Reference  2,  is  more  fre¬ 
quently  and,  in  this  case,  probably  more  appropriately  used.  The 
film  enthalpy  was  used  here,  however,  because  of  its  simplicity 
and  because  it  is  the  value  suggested  for  use  with  the  simplified 
Bartz  equation,  Reference  14.  (Use  of  the  reference  enthalpy  would 
result  in  a  lower  value  of  the  predicted  heat-transfer  coefficient, 
by  about  8  percent.) 
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v*?Jues  of  Prandtl  number  and  viscosity  were  estimated  for  the  two 

~4 

jas  mixtures  and  constant  values  of  Pr  =  0.72  and  u  =  0.5  x  10 
lb/ft- sec  were  used  throughout. 

Retui.«iing  no«  *:?•  *-he  experiments!  results,  Figures  7  and  8, 
the  heat  flux  in  the  tubular  r'^oat  region  increases  with  increasing 
distance  downstream  in  almost  all  cases.  In  most  tests,  a  definite 
time  effect  in  the  value  oi  the  heat-transfer  coefficient  is  apparent. 
This  is  felt  to  be  due  to  the  effect  on  the  calculation  of  enthalpy 
of  the  power  droop  inherent  in  the  battery  power  supply.  The  enthalpy 
calculation  method  currently  being  used  is  such  that  a  somewhat 
lower  than  actual  enthalpy  is  calculated  when  an  appreciable  <;uwer 
droop  occurs  as  is  the  case  at  late  test  times.  it  is  therefore 
felt  that  the  early  time  results,  the  lower  heat-transfer  coeffi¬ 
cient  values,  are  more  representative.  With  the  addition  of  the 
new  diesel  power  supply  this  problem  will  be  eliminated  in  future 
^ests. 

The  experimental  results  are  consistently  below  the  pre¬ 
dicted  heat-transfer  coefficient  in  the  throat  region.  A  line 
which  represents  a  30-percent  decrease  in  the  predicted  value  is 
therefore  included  for  reference.  This  value  is  seen  to  be  in 
good  agreement  with  experiment  (considering,  of  course,  that  the 
early  time  results  are  felt  to  be  more  representative) ,  except 
for  one  test.  Test  No.  948,  Figure  7(a).  No  satisfactory  expla¬ 
nation  for  the  higher  heat-transfer  coefficient  values  in  this 
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test  has  been  found  as  yet.  From  the  results  of  all  other  tests, 
however,  it  is  apparent  that  the  experimental  throat  heat-transfer 
coefficient  is  less  than  the  simplified  Bartz  prediction  by 
approximately  30  percent.  This  has  also  been  observed  in  other 
investigations  of  nozzle  heat  transfer.  The  results  presented 
in  Reference  15  and  also  discussed  in  Reference  2  (see  Fig.  8  of 
that  reference,  particularly)  also  exhibit  experimental  heat- 
transfer  results  which  are  almost  exactly  30  percent  below  the 
simplified  Bartz  prediction  in  the  throat  region. 

In  the  convergent  entry  section  of  the  nozzle,  the  experi¬ 
mental  heat-transfer  coefficient  is  consistently  above  the  predic¬ 
ted  value.  This  reversal  in  the  trend  found  in  the  throat  region 
is  felt  to  be  due  to  radiation  from  the  arc.  This  segment  is 
able  to  "see"  the  arc  and,  therefore,  w.ii^  receive  a  radiant  heat 
input  from  it.  The  radiant  flux  is  apparently  high  enough  to 
cause  a  significant  increase  in  the  apparent  heat-transfer  coeffi¬ 
cient.  No  attempt  to  estimate  this  radiative  flux  and  to  compare 
it  with  the  experimental  results  has  been  made  as  yet.  Unfortu¬ 
nately.  an  adequate  description  of  the  radiation  source,  an  electric 
arc  rotating  at  high  speed  in  an  annular  region,  is  rather  diffi¬ 
cult.  Some  effort  at  estimating  this  radiative  flux  is  expected 
to  be  made  in  the  next  quarter,  however. 

The  experimental  heat-transfer  coefficient  in  the  diver¬ 
gent  exit  section  of  the  nozzle  is  well  below  the  simplified  Bartz 
prediction  in  all  cases.  This  trend  is  as  expected  from  the  com¬ 
parison  in  the  throat  region.  However,  the  experimental  values 
are  consistently  less  than  50  percent  of  the  predicted  values  and 
not  approximately  70  percent  as  found  in  the  throat.  These  even 
lower  coefficients  are  probably  due  to  flow  separation  upstream 
of  the  nozzle  exit.  For  all  three  nozzle  sizes,  the  required 
static  pressure  at  the  exit  for  complete  expansion  is  less  than 
the  actual  static  pressure  of  1  atmosphere.  Flow  separation  would 
therefore  be  expected  upstream  of  the  exit,  with  the  attendant 
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lower  heat  transfer  in  the  separated  region.  Actually,  the  point 
of  separation  should  be  further  upstream  from  the  exit  with  the 
0. 5-inch-diamet jr  throat  nozzle  and  move  downstream  in  going  first 
to  the  0.4  nozzle  and  then  to  the  0.3.  The  discrepancy  with  theory 
should  therefore  become  greater  as  the  throat  diameter  increases, 
and  this  is  indeed  the  case,  as  seen  from  the  figures. b 

3. 2.4.4  Materials  tests 
3. 2. 4. 4.1  Introduction 

The  test  results  for  the  24  graphite,  graphite  phenolic, 
and  silica  phenolic  nozzle  firings  performed  during  this  report 
period  as  a  direct  part  of  the  chemical-erosion  program  are  pre¬ 
sented  below.  Cut-away  photographs  of  the  fired  nozzles  are  pre¬ 
sented  to  reveal  the  nozzle  material  surface  condition  and,  where 
applicable,  char  layer  after  test.  Nozzle  performance  is  presented 
in  terms  of  instantaneous  surface  recession  and  surface-recession 
rate  in  the  throat  region.  Wherever  possible,  comparison  of  the 
results  with  theoretically  predicted  recession  is  made.  In  most 
cases,  the  internal  nozzle  temperature  distribution  is  also  pre¬ 
sented,  and  wall  surface  temperatures  are  indicated  where  possible. 
Again,  appropriate  comparisons  are  made  with  theory. 

Before  turning  to  these  results,  however,  a  brief  dis¬ 
cussion  of  the  test  program  is  appropriate.  The  test  firings  which 
were  outlined  previously  in  Table  II  complete  the  schedule  with 
Mixtures  4  and  5.  Mixture  4  is  made  up  of  the  oxygen  fraction  in 
the  combustion  products  of  a  typical  solid  propellant  that  is  capa¬ 
ble  of  reacting  with  a  nozzle  wall  material,  plus  inert  gases; 

Mixture  5  is  made  up  of  all  inerts.  Both  gases,  of  course,  simu¬ 
late  the  temperature-enthalpy  variation  of  the  typical  solid  pro¬ 
pellant.  The  reader  is  referred  to  References  1  and  2  for  a  complete 

5 In  all  pressure-distribution  calibration  tests,  the  pressure  tap  in 
the  exit  cone  was  always  upstream  of  che  point  where  separation 
would  be  expected.  The  effect  of  separation,  therefore,  was  not 
apparent  in  the  pressure  distribution  results. 
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development  of  the  simulation  requirements  and  tor  the  gas  mixture 
analysis. 

Tests  were  run  over  a  range  of  total  temperature  and 
initial  plenum  pressure  in  order  to  allow  the  investigation  of 
wall  temperature  and  pressure  on  chemical  reaction.  The  emphasis 
is  on  the  high-temperature  condition,  Tq  =  3500*'  K,  however,  since 
this  is  the  approximate  total  temperature  for  the  combustion  prod¬ 
ucts  of  many  solid  propellants. 

A  30-second  firing  duration  was  used  in  almost  all  tests. 
In  the  two  silica  phenolic  tests  at  the  high-temperature  condition, 
this  time  was  shortened  to  approximately  20  seconds  because  of  the 
large  surface  recession.  Also  in  some  tests  premature  arc-out 
occurred  because  of  the  inherent  power  supply  droop  being  unusually 
excessive. 

The  instantaneous  surface  recession  was  calculated  as  a 
function  of  time  for  each  nozzle  firing  from  the  time  variation  of 
plenum  pressure.  This  method  is  discussed  in  detail  in  Reference  2. 
As  indicated  in  Reference  z,  the  expansion  cone  in  the  test  nozzles 
was  felt  to  be  limiting  the  usefulness  of  this  technique  due  to  the 
movement  of  the  minimum  diameter  point,  the  actual  nozzle  throat, 
into  the  expansion  cone  as  surface  recession  progressed.  This 
difficulty  was  eliminated  by  eliminating  the  exit  expansion  cone 
so  that  the  nozzle  exit  was  the  downstream  end  of  the  tubular 
throat.  After  this  fix,  the  agreement  between  measured  and  calcu¬ 
lated  values  of  surface  recession  in  the  tubular  throat  was  excel¬ 
lent,  as  is  apparent  from  the  results  presented  below.  The  new 
test  nozzle  geometry  is  presented  in  Figure  9.  The  internal  nozzle 
wall  profile  both  before  and  after  test  was  determined  by  travers¬ 
ing  the  wall  with  a  pivot -arm,  dial-gauge  micrometer.  Local  surface 
recession,  both  axial  and  circumferential,  was  determined  from 
these  measurements. 

In  most  of  the  firings,  the  internal  wall  temperature 
history  was  measured.  Thermocouples  were  located  in  the  tuoular 
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throat  region  at  three  radial  positions  m  a  plane  normal  to  the 
nozzle  centerline.  The  instrumentation  technique  is  discussed  m 
Reference  2.  The  actual  thermocouple  positions  were  measured  by 
X-ray  photography  as  discussed  in  Section  3  2.3.  Originally,  instru¬ 
mentation  in  three  planes  in  the  tubular  throat  region  was  antici¬ 
pated.  This  was  uoandoned,  however,  since  preliminary  firings 
indicated  little  axial  variation  of  temperature  and  because  of  che 
fairly  jizeable  time  and  cost  involved  in  instrumentating  the 
nozzles.  For  further  details  of  the  overall  testing  technique  and 
other  background  material,  the  reader  is  referred  to  References  1 
and  2. 


3. 2. 4. 4. 2  Post-test  observations 

The  response  of  the  test  nozzles  to  the  simulated  rocket 
environment  was  assessed  in  part  from  post -test  inspection  and 
measurements.  These  results  are  pre?-.'.\\ed  in  this  section.  Photo¬ 
graphs  of  the  fired  nozzles  after  being  sectioned  are  presented  for 
most  tests.  These  photographs  reveal  the  extent  and  uniformity  of 
surface  recession,  the  surface  condition,  and,  in  the  case  of  the 
phenolic3,  the  char  layer.  The  values  of  surface  recession  deter¬ 
mined  from  pre-  and  post-test  measurements  are  also  presented  for 
the  tubular  throat  region.  In  the  case  of  the  phenolics,  the  char 
depths  are  also  discussed. 

The  cut-away  photographs  of  the  test  nozzles  after  firing 
are  shown  in  Figures  10  to  17.  The  original  nozzle  contour  before 
firing  is  indicated  in  each  figure  for  reference.  These  contours 
-re  only  approximate  because  of  alignment  and  enlargement  problems 
in  the  photographic  printing  of  the  composite  of  nozzle  cut-away 
and  original  contour.  Photographs  are  not  presented  for  three  of 
the  fired  nozzles.  In  Test  No.  959  no  surface  recession  occurred 
and,  since  Test  No.  958  was  identical  to  it,  it  was  decided  to 
save  the  nozzle  for  a  future  test.  In  Test  No.  999,  graphite 
phenolic  nozzle  E34,  a  catastrophic  delamination  occurred,  apparently 
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because  of  improper  curing.  Tests  on  the  nozzle  are  now  under  way 
to  check  this  out,  and  therefore  it  is  not  available  for  inclusion 
here.  The  test  conditions  to  which  nozzle  E23  were  exposed.  Test 
Nos.  1017  and  1020,  were  somewhat  erratic,  and  therefore  the  results 
are  questionable  and  are  now  under  study.  Because  of  this,  no  photo¬ 
graph  is  presented  for  this  nozzle. 

The  results  available  from  post-test  observations  and 
measurement  are  discussed  from  the  photographs  and  from  tabulations 
of  the  measured  surface  recessions.  Surface  recession  results  are 
discussed  first,  followed  by  other  observations  and  analysis  of  sur¬ 
face  condition  and  char  depth. 

Absolutely  no  surface  recession  occurs  with  graphite, 
both  Graphitite  GX  and  AT J ,  in  the  inert  gas  environment,  Mixture  5, 
Figures  14  and  15.  Remember  that  this  environment  closely  duplicates 
all  aspects  of  the  solid-propellant  rocket  environment  except  its 
chemical  character.  Also  for  graphite  phenolic.  Figure  16,  the  same 
behavior  is  observed. 

For  Mixture  4,  which  contains  oxygen,  surface  recession 
does  occur  with  graphite  and  graphite  phenolic,  and  its  magnitude  is 
dependent  on  temperature  level,  Figures  10(a)  to  (e)  and  12(a)  to 
(c) .  The  measured  surface  recession  for  graphite  (Graphitite  GX)  is 
summarized  as  follows: 


Graphite,  Graphitite  GX 


Nominal  Total 
Temperature 

(°  K) 

Surface  Recession 
(mils) 

Mixture  4 

Mixture  5 

3500 

74(939) >73(94Q)  >47(997) 

0 (957  &  958) 

2750 

37 ( 1039) 

* 

2000 

0 (1023) 

* 

Numbers  in  parentheses  indicate  test  numbers. 

No  firings  performed  since  no  surface  recession  occurred 
in  the  most  severe  case. 
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All  results  presented  are  for  approximately  the  same  firing  times, 

30  seconds  (see  Table  II) .  The  comparison  between  Mixtures  4  and  5 
indicates  that  the  surface  recession  for  Mixture  4  is  due  solely 
to  chemical  reaction.  The  decrease  in  surface  recession  with 
decreasing  temperature  is  apparently  due  to  kinetic  control  of  the 
oxygen-carbon  reaction  at  the  wall  for  the  two  lower  temperatures. 
This  is  discussed  further  in  Section  3. 2. 4. 4. 4  following.  One  test 
at  the  high-temperature  level  produced  results  which  are  out  of  line 
with  the  others  (see  the  above  table) .  The  test  netup  was  identical 
and  the  test  conditions  were  almost  identical  for  the  three  firings. 
No  obvious  explanation  for  the  discrepancy  is  readily  apparent.  The 
nozzles  could  have  come  from  different  billets  but,  hopefully,  no 
such  variation  between  billets  would  actually  occur. 

As  mentioned  previously,  the  graphite  phenolic  nozzles 
exhibited  no  surface  recession  with  Mixture  5,  the  inert  environ¬ 
ment,  Figure  16.  In  fact,  a  buildup  of  the  solid-resin  pyrolysis 
products  on  the  surface  actually  resulted  in  a  slight  decrease  in 
throat  diameter.  A  significant  weight  loss  occurs,  however,  due 
to  the  resin  decomposition  and  out-gassing.  The  surface  recession 
for  graphite  phenolic  in  the  oxidizing  environment,  Mixture  4,  is 
significant  over  the  entire  temperature  range,  Figure  12.  The 
measured  surface  recession  is  summarized  below: 


Graphite  Phenolic,  MX  4500 


Nominal  Total 
Temperature 

(°K) 

Surface  Recession 
(mils) 

Mixture  4 

h 

Mixture  5 

3500 

64(1001) 

0 ( 1006) 

2750 

54(1040) 

0  ( 1043) 

2000 

4:vxC''4) 

0 ( 1026  &  1028) 

Numbers  in  parentheses  indicate  test  numbers. 

Actually  a  slight  surface  buildup  occurred  in  each  case. 
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Por  Mixture  4,  the  surface  recession  remains  high  with  decreasing 
temperature,  probably  due  to  the  large  number  of  sites  available 
for  chemical  reaction  in  the  "porous"  char  surface.  Again  the  compari¬ 
son  between  Mixtures  4  and  5  points  out  the  significance  of  chemical 
reaction  -  here  it  is  also  apparently  the  only  contributor  to  surface 
recession. 

The  performance  of  silica  phenolic  is  dependent  on  both 
the  temperature  level  and  the  gas  environment,  Figures  13  and  17. 
Surface  recession  is  due  to  "melting"  of  the  silica  and  to  chemical 
reaction  with  the  wall.  A  comparison  of  the  measured  surface 
recession  is  presented  below: 


Silica  Phenolic,  MX2600 


Nominal  Total 
Temperature 

(°K) 

Surface  Recession 
(mils) 

Mixture  4 

Mixture  5 

3500 

160  (94i)  (19  seconds) 

124(960) (21  seconds) 

2750 

60 (1038) 

0  ( 1042) 

2000 

- 

0  ( 1031) 

Numbers  in  parentheses  indicate  test  numbers. 

Note  that  the  recession  at  the  high-temperature  condition  is  for 
firing  tiroes  less  than  30  seconds.  At  this  condition,  the  surface 
recession  is  very  large  -  completely  out  of  line  with  the  perfor¬ 
mance  of  the  other  materials.  Its  degradation  at  high  temperature 
is  largely  due  to  the  "melting"  of  the  silica.  With  Mixture  5, 
the  recession  is  apparently  due  to  melting  only,  whereas  with  Mix¬ 
ture  4  chemical  reaction  is  also  contributing  to  the  surface 
recession.  At  the  intermediate  temperature  condition,  no  melting 
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appears  to  occur  and,  from  a  comparison  of  the  results  for  the  two 
gas  environments,  the  surface  recession  with  Mixture  4  is  apparently 
due  to  chemical  reaction  only. 

The  test  results  presented  and  discussed  above  demonstrate 
conclusively  the  importance  of  chemical  reaction  in  the  erosion  of 
nozzle  wall  materials.  The  surface  recession  for  graphite  and 
graphite  phenolic  was  apparently  due  solely  to  chemical  reaction. 

No  recession  occurred  in  the  inert  environment.  Chemical  reaction 
was  also  a  significant  contributor  to  recession  for  silica  phenolic. 

Post-test  inspection  of  the  nozzle  surface  condition  and 
appearance  was  made  in  each  case  and  revealed  some  interesting 
behavior.  In  the  case  of  the  phenolics  the  surface  condition  and 
appearance  are  dependent  on  both  the  chemical  nature  of  the  gas 
environment  and  the  temperature.  With  Mixture  4,  both  the  graphite 
phenolic  and  silica  phenolic  exhibit  a  "clean"  surface,  apparently 
produced  by  the  reaction  of  the  oxygen  with  the  carbon  produced  in 
resin  pyrolysis.  With  Mixture  5,  however,  the  carbon  products 
remain  in  the  interstices  of  the  char  surface,  apparently  because 
of  the  absence  of  oxygen  in  the  free  stream. 

The  nozzle  wall  for  graphite  phenolic  with  Mixture  4 
becomes  more  rough  and  pitted  as  the  total  temperature  is  decreased, 
Figure  12.  Local  delaminations  are  also  apparent,  particularly  in 
the  tubular  throat  region.  With  Mixture  5,  the  graphite  phenolic 
nozzles  exhibit  a  smooth  wall  surface,  Figure  16.  Some  local 
delaminations  are  also  apparent,  however.  In  all  cases  with  Mix¬ 
ture  4,  the  axial  and  circumferential  variation  of  recession  was 
somewhat  nonuniform. 

For  silica  phenolic  and  Mixture  4  at  high  total  tempera¬ 
ture,  resolidified  silica  or  a  resolidified  silicon  compound  is 
apparent  on  the  surface,  particularly  in  the  subsonic  entry  region. 
Figure  13(a) .  This  is  not  evident  in  the  nozzle  fired  with  Mix¬ 
ture  5,  Figure  17(a).  One  possible  explanation  for  this  is  that 
in  Mixture  5  the  carbon  due  to  resin  pyrolysis  reacts  with  the 
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molten  silica  to  form  carbon  monoxide  and  silicon.  The  silicon, 
due  to  its  low  viscosity,  is  then  removed  mechanically  upon  forma¬ 
tion.  In  Mixture  4,  however,  the  available  carbon  reacts  with  the 
free-stream  oxygen  more  readily  than  with  the  silica  which  apparently 
is  sufficiently  viscous  to  resist  at  least  to  some  extent  the  shear 
forces  which  tend  to  remove  it  from  the  nozzle. 

In  the  graphite  nozzle  firings  in  which  the  surface  reces¬ 
sion  is  appreciable,  the  surface  is  somewhat  rough,  indicating  pre¬ 
ferred  sites  for  chemical  reaction,  for  example,  Figures  10  and  11. 
This  is  particularly  apparent  in  the  ATJ  nozzle.  In  all  cases  where 
surface  recession  occurred,  the  minimum  recession  occurred  at  the 
entry  to  the  tubular  throat. 

The  cut-away  photographs  also  reveal  the  char  layer  in 
the  phenolic  nozzles.  Figures  12,  13,  16,  and  17.  The  char  depths 
are  considerably  greater  for  graphite  phenolic  than  for  silica 
phenolic,  about  a  factor  of  3  for  the  30  second  firing  times.  This 
is  reasonable,  considering  the  lower  thermal  diffusivity  of  the 
silica  phenolic.  The  lower  diffusivity  results  in  a  steeper  tempera¬ 
ture  gradient,  and  therefore  the  resin  pyrolysis  temperature  occurs 
closer  to  the  wall  surface. 

In  all  firings,  the  char  depth  in  the  subsonic  entry 
section  is  the  same  as,  or  greater  than,  that  in  the  tubular  throat 
region.  This  geems  somewhat  anomalous  at  first  sight  since  the 
heat  flux  is  much  lower  in  the  entry  region.  This  behavior  ca..  be 
explained,  however,  by  considering  the  two-dimensional  effect  associ¬ 
ated  with  the  small  diameter  throat  and  the  insensi livity  of  internal 
tenperatuve  to  changes  in  heat  flux  for  low-conductivity  materials. 
The  char  depth  is  also  nearly  independent  of  total  temperature,  for 
example.  Figure  12*  Again,  this  is  due  to  the  low  conductivity  of 
the  phenolic  materials  which  results  in  only  small  internal  tempera¬ 
ture  changes  with  significant  changes  in  surface  heat  flux. 

In  the  silica  phenolic  nozzles,  a  definite  resin  decompo¬ 
sition  region  is  apparent  next  to  the  char  layer.  This  appears  as 
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a  lighter  region  in  the  photographs,  Figures  13  and  17,  and  is 
approximately  half  the  thickness  of  the  char  layer.  No  such 
region  is  apparent  with  the  graphite  phenolic. 

The  char-virgin  plastic  interface  recession  rate  was 
also  assessed  for  the  silica  phenolic  nozzle  firings  from  the 
thermocouple  data  and  the  post-test  measurements  of  the  nozzles. 
These  results  are  discussed  in  Section  3. 2. 4. 4. 5. 

3. 2.4.4. 3  Surface  recession  and  surface-recession  rate  histories 

The  discussions  of  the  preceding  section  were  based 
strictly  on  before-  and  after-the-fact  observations  and  measure¬ 
ments.  During  the  firings  in  which  surface  recession  occurred, 
however,  the  surface  recession  rate  varied  over  wide  ranges.  A 
significant  transient  time  exists  early  in  each  firing  before  a 
quasi-steady  recession  condition  is  achieved.  Also,  once  the 
quasi-steady  condition  is  reached,  the  surface-recession  rate 
remains  a  variable,  decreasing  with  time  due  to  the  decreasing 
heat-  and  mass-transfer  coefficients  associated  with  the  increas¬ 
ing  throat  diameter.  These  considerations  are  particularly 
important  in  making  general  conclusions  regarding  surface  reces¬ 
sion  and  surface-recession  rate.  An  average  surface-recession 
rate  calculated  only  from  pre-  and  post-test  measurements  and 
the  firing  duration  can  be  quite  misleading  if  it  is  applied  for 
firing  times  or  geometries  other  than  those  for  which  the  results 
were  obtained.  A  knowledge  of  the  instantaneous  surface  recession 
and  surface-recession  rate  is  therefore  important.  These  results 
are  presented  in  this  section. 

The  surface  recession  and  surface-recession  rate  are 
plotted  against  time  in  Figures  18  to  22  for  all  nozzle  firings 
in  which  surface  recession  occurred  at  the  high-temperature  con- 
dition.  The  data  for  the  firings  at  the  intermediate  and  low 

6As  discussed  briefly  in  Section  3.2.3  and  in  detail  in  Reference  2, 
the  instantaneous  surface  recession  and  surface-recession  rate 
are  calculated  from  the  measured  plenum  pressure  decay. 
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temperature  conditions  ware  in  the  reduction  process  at  the  tine 
of  writing  and  therefore  are  not  presented.  All  plots  are  to  the 
same  scale  for  ready  comparison.  Tho  test  conditions,  including, 
where  available,  the  wall  surface  temperature  indicated  by  the 
pyrometer,  and  a  comparison  of  the  measured  and  calculated  final 
throat  diameter,  are  presented  in  each  figure.  The  agreement 
between  these  throat  diameters  is  good  in  every  case.  The  results 
presented  are  therefore  felt  to  represent  closely  the  actual 
instantaneous  variation  of  surface  recession  and  surface-recession 
rate  in  the  tubular  throat  region  during  each  firing. 

The  theoretically  predicted  surface-recession  rate  ie 
also  presented  for  the  graphite-nozzle  firings  with  Mixture  4. 

The  prediction  technique  and  the  comparison  of  the  experimental 
results  with  prediction  are  discussed  in  Section  3. 2. 4. 4. 4  follow¬ 
ing. 

The  surface  recession  exhibits  two  different  types  of 
variation  with  time,  depending  on  the  time  during  the  firing.  Fig¬ 
ures  18  to  22. .  At  early  tildes,  both  the  surface  recession  and  the 
surface-recession  rate  are  increasing.  At  later  times  the  surface 
recession  is  of  course  still  increasing  but  the  surfiice-recession 
rate  is  decreasing.  The  early  time  variation  is  a  transient 
response  to  the  step  change  in  conditions  at  the  start  of  firing. 
After  the  nozzle  has  accommodated  itself  to  this  step  change,  the 
second,  later-time  variation  is  apparent.  The  surface  recesrion 
in  this  regime  may  be  termed  quasi-steady  in  that,  even  though  the 
surface-recession  rate  Is  decreasing,  this  variation  occurs  only 
after  the  equilibration  to  the  step  change  in  conditions  at  the 
start  of  firing  has  occurred.  Note  that  as  much  as  18  seconds 
may  be  required  to  reach  the  quasi-steady  condition  for  the  range 
of  conditions  end  materials  presented. 

The  surface-recession  rates  for  all  graphite-nozzle 
firings  at  tha  high-temperature  condition  are  presented  in  Fig¬ 
ure  23.  The  three  Graphitite  GX  graphite  firings  at  high  pressure, 
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Test  Nos.  939,  940,  and  997,  show  reasonable  agreement  and,  as  indi¬ 
cated  in  Figure  24,  fall  within  a  fairly  narrow  band  with  a  maximum 
surface-recession  rate  at  early  time  of  about  4  mils/sec  and  a  mini¬ 
mum  at  the  end  of  the  firings  of  something  higher  than  1.5  mils/sec. 
In  Figure  23,  the  ATJ  graphite  firing.  Test  No.  1000,  is  seen  to 
have  a  quite  different  trend  than  the  three  Graphitite  GX  firings 
at  the  same  nominal  conditions.  The  surface-recession  rate  is 
lower  at  early  times  but  levels  off  such  that  at  the  end  of  the 
firings  it  has  a  recession  rate  comparable  to  that  of  the  three 
Graphitite  GX  firings.  This  early-time  behavior  is  probably  due, 
at  least  in  part,  to  the  higher  thermal  diffusivity  of  ATJ  graphite 
at  low  temperature  which,  in  turn,  results  in  a  lower  wall  surface 
temperature  at  early  times.  If  the  firings  had  not  been  terminated 
at  30  seconds,  it  appears  that  the  recession  rate  for  ATJ  would  be 
higher  than  that  for  Graphitite  GX  at  times  in  excess  of  30  seconds. 

By  decreasing  the  initial  plenum  pressure,  (Test  No.  998, 
Fig.  23) ,  and  therefore  the  heat-  and  mass-transfer  coefficients  and 
probably  the  wall  temperature,  a  significant  decrease  in  surface- 
recession  rate  occurs  -  better  than  a  factor  of  2  for  a  reduction 
in  initial  plenum  pressure  from  300  psia  to  177  psia.  The  rate  is 
also  almost  invariant  with  time,  at  least  for  the  30  seconds  of 
firing  duration.  A  further  discussion  of  the  comparison  between 
these  results  is  presented  in  the  discussion  of  the  predicted 
surface-recession  rates,  Section  3. 2. 4. 4. 4. 

A  comparison  of  the  recession  rates  for  graphite,  graphite 
phenolic,  and  silic.a  phenolic  is  presented  in  Figure  24  for  the 
high-temperature  and  high-pressure  condition.  The  graphite  tests 
are  presented  by  the  cross-hatched  band  which  encloses  the  results 
of  the  three  Graphitite  GX  firings.  The  surface-recession  rate 
for  graphite  phenolic  falls  within  the  band  for  the  graphite  firings 
except  at  late  ti:r>es  when  it  is  lower.  At  early  times  it  is  close 
to  the  upper  lirr.it  for  graphite.  The  silica  phenolic  recession-  rate 
is  much  higher  than  that  for  graphite  and  graphite  phenolic.  The 
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maximum  value  is  about  12.5  mils/sec  after  5  seconds  of  firing, 
the  start  of  quaai-steady  conditions.  Note  that  a  quasi-steady 
recession  is  achieved  much  earlier  in  the  firing  than  for  the 
other  materials,  because  of  its  low  thermal  diffusivity.  At 
If*  seconds,  ths  recession  rate  is  about  twice  that  of  the  other 
materials,  7.5  mils/sec.  In  this  firing,  the  throat  diameter 
increased  by  a  factor  of  better  than  2,  and  therefore  the  area 
by  better  than  4,  over  the  20-second  firing  duration. 

The  results  of  an  early  firing  on  an  ATJ  graphite 
nozzle.  Test  No.  795,  were  presented  in  Reference  2.  This  firing 
was  performed  before  the  swirl  component  was  eliminated  from  the 
arc-plasma-generator  flow,  however.  The  surface-recession  rate 
for  this  firing  is  presented  in  Figure  25,  together  wit';  the 
results  of  the  recent  ATJ  firing,  Test  No.  1000.  The  early  test 
before  swirl  was  eliminated  has  a  substantially  higher  surface- 
recession  rats.  This  is  felt  to  be  due  to  the  higher  heat-  and 
mass-transfer  coefficients  caused  by  swirl  in  the  early  test. 

The  higher  heat-transfer  coefficient  can  be  expected  to  result 
in  a  higher  wall  temperature  at  an  earlier  time  in  the  firing, 
thus  contributing  to  a  higher  surface-recession  rate  when  the 
recession  is  reaction  rate  controlled.  The  higher  mass-transfer 
coefficient  results  in  more  oxygen  available  for  reaction  at  the 
surface  and,  therefore,  a  higher  surface-recession  rate  when 
recession  is  diffusion  rate  controlled,  that  is,  when  equilibrium 
obtains.  The  slightly  more  severe  chamber  conditions  in  the 
earlier  firing  is  a  further  contribution  to  the  discrepancy  in 
the  results.  The  primary  contributor,  however,  is  felt  to  be 
the  swirl. 


3. 2. 4. 4. 4  Comparison  of  experiment  with  theoretical  prediction 

Before  discussing  the  comparison  of  experimental  sur¬ 
face  recession  rate  with  prediction,  a  brief  discussion  of  the 
meaning  of  such  a  comparison  is  in  order.  The  predicted  surface 
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recession  rate  is  based  on  the  assumption  that  all  the  reactant, 
in  this  case  oxygen,  that  diffuses  to  the  wall  reacts  with  it  (no 
mechanical  erosion  is  considered) ,  This  rate  of  diffusion  of  the 
reactant  is  given  by  the  term 
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At  low  enough  wall  temperatures,  however,  the  surface  recession 
is  influenced  by  the  reaction  rate  as  well  as  the  diffusion  rate. 

In  this  case,  not  all  the  reactant  available  does,  in  fact,  react. 
The  experimental  surface  recession  will  then  be  lower  than  that 
predicted.  The  reaction  rate  is  generally  expressible  by  an 
equation  of  the  Arrhenius  form  ^jPj  =  AjPj  e^J^ ,  where  the 
units  are  the  same  as  those  for  the  diffusion  rate  and  where  Pj 
and  T  are  the  partial  pressure  of  reactant  and  the  temperature. 
Note  that  the  reaction-rate  constant  is  particularly  sensitive  to 
temperature  in  that  it  appears  as  an  exponential.  In  summary, 
the  surface  recession  corresponding  to  diffusion  rate  control 
represents  the  upper  bound  on  chemical  erosion.  A  surface 
recession  lower  than  this  value  indicates  that  the  wall  reaction 
is  reaction  rate  (kinetically)  controlled. 

An  experimental  surface-recession  rate  higher  than 
the  equilibrium  prediction  indicates  mechanical  erosion,  since 
our  prediction  is  the  upper  limit  on  chemical  reaction.  Mechanical 
erosion  actually  can  present  an  insidious  complication  in  com¬ 
paring  experiment  with  prediction  since  any  additional  surface 
recession  due  to  mechanical  erosion  tends  to  compensate  for  any 
decrease  due  to  kinetic  control  of  the  surface  reaction.  Fortu¬ 
nately,  no  mechanical  erosion  of  graphite  or  graphite  phenolic 
occurred  with  Mixture  5  and,  therefore,  would  not  be  expected  in 
the  chemically  reactive  mixtures,  unless,  however,  the  chemical 
reaction  weakens  the  wall  material  and  allows  mechanical  erosion 
to  take  place.  An  appropriate  test  program  could  of  course. 
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determine  the  presence  or  absence  of  mechanical  erosion  in  the 
chemical  environment,  but  is  not  included  as  part  of  the  chemical- 
erosion  program. 

The  derivation  of  the  surface-recession  rate  prediction 
was  presented  in  Reference  2  and  is  not  repeated  here.  The  mass- 
transfer  coefficient  required  for  the  calculation  was  determined 
from  the  simplified  Bartz  equation  for  heat-transfer  coefficient 
and  the  analogy  between  heat  and  mass  transfer 

p'  ueCM'  =  P  '  uecH '  ^,2/3 

where  Le  is  the  Lewis  number.  Unity  was  chosen  as  a  convenient 
yet  reasonable  value  for  the  Lewis  number.  From  the  heat-transfer 
calibration  tests,  it  was  found  that  the  experimental  heat-transfer 
coefficient  was  approximately  30  percent  lower  than  the  Bartz  pre¬ 
diction.  The  experimental  mass-transfer  coefficient  should  there¬ 
fore  also  be  lower  than  the  predicted  value  by  the  same  amount. 

Two  theoretical  lines  are  indicated  on  the  plots,  one  based  on  the 
Bartz  predicted  mass-transfer  coefficient  and  the  other  30  percent 
lower,  which  is  based  on  the  heat-transfer  calibration  tests.  The 
lower  curve  should  be  closer  to  fact  for  chemical  equilibrium  with 
the  wall. 

The  experimental  surface-recession  rate  is  compared  with 
the  predicted  value  for  all  graphite-nozzle  firings  in  Figure  18 
and  is  presented  in  composite  form  for  the  three  Graphitite  GX 
firings  at  the  same  nominal  conditions  in  Figure  26.  At  these 
conditions,  the  experimental  recession  rate  is  slightly  lower 
than  the  predicted  value,  departing  by  the  largest  amount  at  late 
times  when  the  pressure  and  possibly  the  wall  temperature  have 
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decreased  due  to  the  enlargement  of  the  throat.7  At  early  times, 
therefore,  the  surface  recession  is  very  nearly,  if  not  completely, 
controlled  by  the  diffusion  of  oxygen  to  the  wall.  At  the  late 
times,  however,  the  wall  temperature  apparently  dropped  to  a  suffi- 
cently  low  level  that  the  reaction  at  the  wall  is  reaction  rate 
influenced. 

At  the  high-temperature  and  low-pressure  condition,  Fig¬ 
ure  18(d),  the  experimental  surface-recession  rate  is  less  than 
half  the  predicted  rate.  The  final  throat  diameter  (and  mass-flow 
rate)  and,  therefore,  the  final  free-stream  conditions  in  this 
firing  are  almost  identical  to  that  at  the  end  of  firing  in  the 
tests  of  Figures  18(a)  to  (c) .  The  difference  in  surface-recession 
rate  and  the  comparison  with  theory  between  these  two  conditions 
must  therefore  be  due  to  a  difference  in  wall  temperature.  Because 
of  the  lower  heat-transfer  coefficient  in  the  firing  of  Figure  18(d) , 
the  wall  temperature  would  not  be  expected  to  be  as  high  as  that 
for  the  firings  of  Figures  18(a)  to  (c) .  A  lower  reaction  rate 
and  therefore  a  lower  surface-recession  rate  obtains.  The  compari¬ 
son  of  the  two  cases  will  be  investigated  further  with  the  one- 
moving  -boundary  computer  program. 

Kinetic  or  reaction  rate  influence  is  also  evident  in 
the  two  graphite  firings  at  lower  total  temperature  for  which  the 
data  reductions  are  not  complete.  In  the  lowest  total  temperature 


7Aithough  a  decrease  in  wall  temperature  at  later  times  in  a  firing 
seems  a  bit  out  of  order  at  first  glance,  it  is  indeed  a  real 
possibility.  With  a  sufficiently  large  decrease  in  heat-transfer 
coefficient,  a  decrease  in  wall  temperature  can  occur  in  the 
absence  of  surface  recession,  even  though  the  wall  recovery 
enthalpy  remains  constant.  The  one-moving-boundary  computer  pro¬ 
gram  will  be  used  to  checkout  this  possibility  for  the  case  of 
chemical  erosion.  Note  that,  as  seen  in  the  plots  of  temperature 
history,  Figure  27,  the  wall  temperature  indicated  by  the  pyrometer 
does  decrease  with  firing  time  at  late  times. 

0 

Note  that  the  pyrometer  indication  of  wall  temperature  does  not 
indicate  a  substantial  difference  between  the  two  conditions. 
However,  as  indicated  previously,  there  is  some  doubt  in  the 
validity  of  the  pyrometer  indications  due  to  th^.  effects  of  looking 
through  the  luminous  exit  flow. 
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firing,  nominally  2000°  K,  the  experimental  surface-recession  rate 
is  zero,  whereas  the  predicted  value  is  a  constant  at  approximately 
6.5  mil3/sec.  The  wall  temperature  in  this  case  will  of  course  be 
substantially  less  than  2000°  K. 

In  the  ATJ  graphite  firing.  Figure  19,  the  experimental 
surface-recession  rate  is  somewhat  below  the  predicted  rate,  but 
approaches  it  at  late  times.  This  low-recession  rate  is  somewhat 
surprising  when  one  considers  the  lower  density  and,  therefore 
the  expected  availability  of  additional  sites  for  chemical  reac¬ 
tion.  However,  it  may  be  due  to  the  higher  thermal  diffusivity 
of  ATJ  graphite  at  low  temperature;  a  low  wall  temperature9  or 
possibly  a  longer  time  to  achieve  a  high  wall  temperature  may 
result.  (Again,  this  will  be  investigated  further  with  the  one- 
moving-boundary  computer  program.)  Due  to  the  apparently  low 
temperature,  the  surface  recession  is  reaction  rate  influenced, 
at  least  up  to  the  30-second  test  time  of  this  firing. 

In  summary,  reaction-rate  control  of  surface  recession 
is  evident,  at  least  to  some  extent,  in  all  of  the  graphite-nozzle 
firings.  The  magnitude  of  this  effect  is  apparently  associated 
primarily  with  the  wall  surface  temperature  rather  than  the  pres¬ 
sure.  At  the  high  wall  temperatures  experienced  at  early  firing 
times  in  the  high  total  temperature,  high  pressure  firings,  the 
surface  recession  is  close  to  being,  if  not  completely,  diffusion 
rate  controlled. 

Further  analysis  of  the  graphite  firing  results  will  be 
performed  with  the  one-moving-boundary  computer  program.  Sec¬ 
tion  2.2,  to  investigate  further  the  observed  experimental  behavior 
of  the  surface  recession  and  to  verify  the  explanations  presented 
above.  The  charring-ablator  computer  program  will  also  be  used 
for  an  analysis  of  the  experimental  results  with  graphite  phenolic 
and,  conversely,  the  experimental  results  will  be  used  as  a  guide 


9Note  that  in  this  case,  the  pyrometer  reading  does  support  this 
explanation  of  the  behavior. 
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in  determining  reaction-rate  data  for  use  in  the  program.  These 
further  analyses  will  also  include  an  investigation  of  the  experi¬ 
mentally  determined  internal  temperature  distributions  discussed 
below. 


3. 2.4.4. 5  Temperature  histories 

The  internal  wall  temperature  histories  and  the  internal 
temperature  distributions  obtained  therefrom  provide  a  further 
means  of  assessing  the  chemical  behavior  at  the  nozzle  wall.  This 
information  can  be  used  in  the  determination  of  wall  surface  tempera¬ 
ture,  which  of  course  affects  chemical  reaction,  and  in  comparisons 
with  the  theoretical  prediction  techniques  being  developed  in  the 
program.  The  temperature  data  required  to  provide  this  information 
were  obtained  by  thermocouple  instrumentation  at  three  radial  loca¬ 
tions  in  a  plane  perpendicular  to  the  nozzle  centerline  that  also 
goes  through  the  center  of  the  tubular  throat  (see  Ref.  2)  .  Three 
radial- thermocouple  locations  were  used  so  that  the  internal  tem¬ 
perature  distribution  at  any  time  in  the  firing  is  available. 

The  interpretation  of  the  results  is  complicated  by  sur¬ 
face  recession,  the  distance  between  the  thermocouples  and  heated 
surface  continually  decreasing,  decreasing  wall  heat  flux  with 
time,  and,  in  the  case  of  the  phenolics,  the  existence  of  a  char 
layer  and  pyrolysis  region  between  the  heated  surface  and  the 
thermocouples.  Because  of  these  problems  it  is  difficult  in  most 
cases  to  obtain  definitive  information  on  surface  temperature  and 
chemical  reaction  directly.  The  primary  utility  of  the  experi¬ 
mental  distributions  is  in  comparison  of  the  experimental  results 
with  those  obtained  by  the  theoretical  prediction  techniques  being 
developed  under  the  contract.  From  these  comparisons,  together  of 
course  with  all  other  available  experimental  information,  the 
appropriate  chemical-reaction  rate  data  for  use  in  the  program  can 
be  determined.  Since  the  one-moving -boundary  computer  program 
(Section  2.2)  has  only  recently  become  available  and  the  charring 
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ablater  computer  program  (Section  2.3)  is  not  yet  complete,  the 
available  temperature  history  results  have  not  yet  been  put  to 
full  use.  The  pertinent  information  now  available  from  the  tem¬ 
perature  histories  is  presented  and  discussed  below. 

The  internal  wall  temperature  histories  for  all  nozzle 
materials  and  both  test  gas  mixtures  are  presented  in  Figures  27  to 
32,  at  the  high-temperature  test  condition.  The  results  for  the 
lower  temperature  conditions  are  in  the  data-reduction  process 
and,  therefore,  are  not  presented.  The  test  conditions  and  the 
locations  of  the  thermocouples  are  indicated  in  each  figure. 10 
The  pyrometer  indication  of  the  wall  temperature  variation  is 
also  indicated  where  available.  For  the  graphite-nozzle  firings 
with  Mixture  5  (no  surface  recession  and  no  chemical  reaction) , 
the  experimental  results  are  compared  with  the  theoretical  tem¬ 
perature  histories  calculated  with  the  computer  program  of  Ref¬ 
erence  16.  This  program  solves  the  differential  equation  for 
the  transient  heating  of  an  axi symmetric  body  by  a  finite  dif¬ 
ference  technique.  The  program  allows  variable  material  thermal 
properties  and  time-dependent  heat-transfer  coefficient  and 
recovery  enthalpy.  The  heat-transfer  coefficient  used  was  that 
predicted  by  the  Bartz  equation  but  reduced  by  30  percent,  as 
is  indicated  appropriate  by  the  heat-transfer  calibration  tests. 

In  the  comparison  of  the  experimental  and  theoretical 
temperature  histories  in  the  nonablating,  nonreacting  cases. 

Figure  30,  agreement  in  both  magnitude  and  trend  for  the  complete 
firing  in  both  tests  is  surprisingly  good.  The  discrepancy 
between  experiment  and  theory  never  exceeds  90°  K.  This  occurs 
at  high  temperature,  and  part  of  it  is  due  to  the  inadequacy  of 
the  linear  temperature-millivolt  thermocouple  calibration  assumed 
for  data  reduction.  From  this  agreement  with  theory,  it  appears 

l0Thermocouple  locations  in  parentheses  indicate  nominal  values; 
an  X-ray  photograph  of  this  nozzle  was  not  taken  due  to  an 
oversight. 
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that  the  thermocouple  instrumentation  technique  does  indeed  yield 
reliable  results. 

The  theoretical  surface  temperatures  are  also  presented 
in  the  figures  and,  based  on  the  agreement  at  the  internal  points, 
they  should  be  accurate.  Unfortunately,  no  pyrometer  data  are 
available  for  comparison  purposes  in  these  tests.  The  pyrometer 
was  set  up  but  no  reading  was  obtained.  Based  on  the  theoretical 
surface  temperature,  no  reading  is  exactly  what  would  be  expected 
-  the  pyrometer  was  set  on  a  high-temperature  scale  that  had  a 
minimum  temperature  above  the  maximum  value  indicated  by  the  theo¬ 
retical  prediction  (this  scale  setting  was  chosen  because  it  was 
the  one  found  to  be  correct  for  the  same  conditions  with  Mixture  4) . 

In  the  comparable  firing  with  Mixture  4  where  both 
surface  recession  and  chemical  reaction  occurred.  Figure  27(a), 
the  temperature  at  the  nearest-to-the  surface  location  rose  con¬ 
siderably  faster  than  that  at  the  same  location  with  Mixture  5. 

This  response  is  felt  to  be  due  both  to  the  chemical  reaction 
occurring  at  the  wall  and  to  the  surface  recession  associated  with 
it.  Thermocouple  1  in  the  firing  of  Figure  27(b)  did  not  appear 
to  be  functioning  properly  and,  therefore,  should  not  be  used  in 
interpreting  the  results. 

In  both  of  these  firings.  Figure  27,  the  pyrometer 
indicates  a  decreasing  surface  temperature  in  the  time  ranges  for 
which  information  is  available.  If  the  pyrometer  is  indeed  indi¬ 
cating  the  actual  variation,  it  may  well  be  due  to  the  decreasing 
heat-  and  mass-transfer  coefficient  with  surface  recession.  This 
trend  will  be  checked  further  with  the  one-moving -boundary  com¬ 
puter  program. 

The  graphite  phenolic  firings  with  Mixtures  4  and  5, 
Figures  28  and  31,  appear  similar  in  the  time  ranges  for  which 
comparable  results  are  available.  Unfortunately,  one  thermocouple 
failed  in  each  firing.  With  Mixture  4,  the  thermocouple  closest 
to  the  surface  reached  a  temperature  quite  close  to  the  surface 


-70- 


temperature  at  indicatad  by  the  pyrometer;  at  the  end  of  test. 
Thermocouple  1  it  only  about  0.040  inch  from  the  surface. 

The  motion  of  the  char  -virgin  plastic  interface  is  not 
apparent  in  any  of  the  therm^cox’ple  traces  in  either  of  the  graphite 
phenolic  firings.  Mo  abrupt  changes  in  temperature  occurred  to 
signal  the  passing  of  the  interface.  However,  both  Thermocouples  1 
and  2  were  in  the  char  layer  at  the  end  of  the  firings. 

The  silica  phenolic-nozzle  firings.  Figures  29  and  32, 
are  particularly  interesting,  at  least  academically.  In  these 
firings,  the  motion  of  the  char-virgin  plastic  interface  is  identi¬ 
fiable  from  the  thermocouple  histories  and  Thermocouple  1  was  over¬ 
taken  by  surface  recession.  From  Figures  29  and  32,  the  variation 
of  temperature  with  time  is  characterized  by  a  slow  increase  in 
temperature  followed  by  an  abrupt  rise  to  high  temperature. lx  This 
is  particularly  evident  in  Figure  29  for  bcch  Thermocouples  1  and  2. 
This  sharp  increase  in  temperature  occurs  within  the  first  1/8-inch 
of  the  nozzle  wall.  The  start  of  the  abrupt  rise  probably  signals 
the  end  of  resin  pyrolysis  and  the  start  of  the  motion  of  the  char 
surface  across  the  thermocouple  location.  Post-test  measurements 
of  the  extent  of  the  char  layer  support  this  interpretation.  Pyrol¬ 
ysis,  therefore,  is  seen  to  occur  at  temperatures  less  than  800°  K 
(1000°  F). 

In  both  firings,  the  surface  recession  overtook  Thermo¬ 
couple  1.  The  point  in  time  at  which  this  occurred  is  indicated 
in  the  figures  and  was  determined  from  the  calculated  surface 
recession  and  the  thermocouple  locations.  The  probable  tempera¬ 
ture  history  up  to  this  point  is  indicated  in  Figure  29. 

The  motion  of  the  char-virgin  plastic  interface  as  deter¬ 
mined  from  the  temperature  histories  is  presented  in  Figures  33 
and  34.  The  plot  on  the  left  is  the  recession  referenced  to  the 

1XXn  Figure  32,  Thermocouple  2  started  an  abrupt  rise  in  tempera¬ 
ture  about  1/2  second  after  the  firing  was  terminated,  followed 
almost  immediately  by  a  drop  in  temperature  as  cool-down  started. 
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location  of  the  wall  surface  at  the  start  of  the  firing.  The  right- 
hand  plot  shows  the  variation  of  char  depth  determined  from  the 
surface  recession  and  the  recession  of  the  char-virgin  plastic 
interface.  The  recession  determined  from  measurement  of  the  inter¬ 
face  location  in  the  sliced  nozzles  after  firing  is  also  indicated. 

The  agreement  between  the  value  determined  from  post-test  measure¬ 
ment  and  that  predicted  from  the  temperature  history  information  is 
good.  The  variation  of  char  depth  with  time  is  characterized  by  a 
rise  to  a  value  which  is  invariant  with  time  in  the  latter  part  of 
the  firing.  This  depth  for  the  firing  with  Mixture  4  is  approximately 
110  mils  and  with  Mixture  5  about  140  mils.  This  compares  with  the 
measured  char  depths  of  110  mils  and  120  mils,  respectively. 

3.2.5  Summary 

The  results  of  a  number  of  nozzle  firings  of  typical  solid- 
propel.lant  rocket-motor  exposed  materials  were  presented  and  dis¬ 
cussed,  together  with  the  results  of  appropriate  heat-transfer  and 
pressure -distribution  calibration  tests.  These  tests  were  per¬ 
formed  as  part  of  the  chemical-erosion  program  and  represent  almost 
half  of  the  testing  to  be  performed  under  this  phase.  The  nozzle 
materials  were  Graphitite  GX  graphite,  ATJ  graphite,  MX4500  graph¬ 
ite  phenolic,  and  MX2600  silica  phenolic.  Tests  were  performed 
with  Mixtures  4  and  5.  both  mixtures  simulating  all  aspects  of  the 
rocket  environment  except  its  chemical  nature.  Mixture  5  is  com¬ 
pletely  inert  chemically,  whereas  Mixture  4  contains  an  amount  of 
oxygen  equivalent  to  that  in  typical  solid-propellant  combustion 
products  that  is  able  to  react.  The  nominal  range  of  test  con¬ 
ditions  was  a  gas  total  temperature  from  3500°  K  to  2000°  K  and  a 
plenum  pressure  from  300  psia  on  down. 

The  nozzle-firing  results  were  presented  in  terms  of  the 
measured  total  surface  recession,  the  time  variation  of  surface 
recession  and  surface-recession  rate,  and  the  internal  temperature 
histories.  Cut-away  photographs  of  most  nozzles  were  also  pre¬ 
sented  and  the  pertinent  information  available  from  them  was 
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discussed.  This  includsd  surface  condition,  uniformity  of  surface 
recession,  and  the  char  layer.  Wherever  possible,  the  experimental 
results  were  compared  with  theory.  The  calibration  test  results 
were  presented  in  the  form  of  the  heat-transfer  coefficient  axial 
distribution  and  the  pressure-ratio  axial  distribution.  Again, 
comparisons  were  made  with  theory. 

A  summary  of  the  results  and  conclusions  is  presented  below: 

(1)  Chenicsl  reaction  appeared  to  be  the  only  contributor 
to  surface  recession  for  graphite  and  graphite  phenolic;  no  sur¬ 
face  recession  occurred  with  these  materials  in  the  inert  environ¬ 
ment,  Mixture  5,  whereas  the  surface  recession  was  significant 
with  Mixture  4,  the  chemically  reactive  environment. 

(2)  Prom  the  comparison  of  experiment  and  theory,  the 
influence  of  both  species  diffusion  and  chemical  kinetics  on  sur¬ 
face  recession  was  apparent  for  the  range  of  test  conditions. 

(3)  At  high  total  temperature,  the  performance  of  graphite 
and  graphite  phenolic  are  comparable;  at  lower  temperatures  graph¬ 
ite  was  superior,  however. 

(4)  Silica  phenolic  exhibited  severe  surface  recession  at 
high  total  temperature  due  to  "melting;"  chemical  reaction  was 
also  a  contributor  in  the  chemical ly-reactive  environment,  Mix¬ 
ture  4. 

(5)  Where  comparison  was  made,  the  experimental  internal- 
temperature  histories  agreed  favorably  with  the  theoretical  pre¬ 
diction. 

(6)  The  measured  heat-transfer  coefficient  was  approxi¬ 
mately  30  percent  below  the  simplified  Bartz  prediction  and  in 
agreement  with  results  of  other  similar  tests* 

The  analysis  of  the  results  is  incomplete  in  that  not  all 
test  data  have  been  reduced  at  the  time  of  writing  and  not  all 
theoretical  comparison  were  able  to  be  made.  This  will  be  per¬ 
formed  during  the  next  quarterly  period. 
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3.3  Particle-Impact  Test  Program 
3.3.1  Introduction 

The  particle-impact  program  is  an  experimental  investi¬ 
gation  of  the  effects  of  metal-oxide  particle  impingement  on  the 
behavior  of  rocket-nozzle  materials.  The  wall  response  and  par¬ 
ticle  behavior  associated  with  particle  impact  is  assessed  in 
terms  of  the  individual  mechanisms  associated  with  or  influenced 
by  particle  impingement.  These  mechanisms  include: 

kinetic  energy  transfer 
internal  thermal  energy  transfer 
chemical  reaction 
shear  due  to  impact 
radiative  heat  transfer 

boundary- layer  perturbations  which  affect: 
convective  heat  transfer 
mass  transfer 
shear 

The  importance  of  these  mechanisms  and  their  contribution  to 
nozzle-material  wall  degradation  will  be  studied  by  two  compli¬ 
mentary  and  concurrent  approaches.  These  are  the  definitive 
measurements  of  both  the  material  response  and  the  particle 
variables  and  behavior  before,  during,  and  after  impact,  and  the 
observation  and  measurement  of  material  response  in  pairs  of  tests 
that  are  identical  except  for  the  absence  of  particles  in  one  and 
their  presence  in  the  other. 

Testing  in  the  particle-impact  program  will  be  started 
early  in  the  next  quarterly  period,  and  before  the  testing  in  the 
chemical-erosion  program  (on  Mixtures  1-3)  is  completed.  These 
tests  will  take  the  form  of  checkout  runs  to  obtain  a  general 
feel  for  the  conduct  of  the  particle-impact  tests  and  to  define 
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any  problem  areas  which  will  require  further  effort  before  launch¬ 
ing  into  the  primary  testing  phase  of  the  program.  If  any  such 
problems  do  come  up,  the  chemical-erosion  program  will  be  picked 
up  again  while  the  appropriate  fixes  are  being  made.  The  effort 
in  the  program  during  the  third  quarter  consisted  of  the  defini¬ 
tion  and  investigation  of  instrumentation,  particularly  for  mea¬ 
surement  of  particle  velocity,  and  the  start  of  the  test  setup 
for  the  program. 

3.3.2  Testing  technique 

Testing  will  be  conducted  under  controlled  laboratory  con¬ 
ditions  which  closely  duplicate  the  solid-propellant  rocket  com¬ 
bustion  products  environment  including  the  presence  of  solid 
particles  or  liquid  droplets  of  aluminum  oxide.  The  Vidya  arc- 
plasma  generator  provides  the  combustion  products  environment. 

The  test  sections  are  typical  rocket -nozzle  wall  materials, 
graphite,  graphite  phenolic,  silica  phenolic,  and  tungsten.  The 
exit  flow  from  the  arc-plasma  generator  impinges  on  the  test 
sections  which  are  placed  downstream  of  the  plasma-generator  exit 
plane  at  any  desired  preset  angle  of  attack.  Briefly,  this  out¬ 
lines  the  test  technique  and  procedure;  the  details  are  presented 
below. 

The  description  of  the  arc-plasma  generator  and  its  oper¬ 
ating  characteristics  is  presented  in  References  1  and  2  and  is 
not  discussed  further  here.  The  reader  is  referred  to  these 
references  for  this  background  information.  The  test  gas  to  be 
used  in  the  program  for  most  tests  is  Mixture  4. 12  A  number  of 
tests  are  also  anticipated  with  Mixture  5,  which  consists  of 
inert  gases  only,  in  order  to  investigate  particle  effects  in 
the  absence  of  chemical  reaction  of  the  wall  material  with  the 


18Again,  the  reader  is  referred  to  References  1  and  2  for  the 
background  information  on  the  test  gas  mixtures. 
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test  gas.  The  particles  are  high  purity  alumina,  Al^O^,  Norton 
Company  electronic  grade  Alundum,  99.65  percent  Al203  minimum. 
Two  different  particle  sizes  will  be  used  -  one  for  which  the 
average  size  on  a  weight  basis  is  12  microns,  and  the  other 
5  microns.  The  size  distributions  provided  by  the  manufacturer 
are  indicated  in  the  following  sketch: 


Z  0  I _ 1— — I  I  1  1  1  1 .11 

fcj  1  2  5  10 

Diameter,  microns 


The  distributions  and  individual  particle  shapes  will  also  be 
checked  by  electron  microscope.  Since  the  test  procedure  results 
in  all  particles  in  the  exit  flow  impacting  on  the  test  sample, 
mass  loadings  typical  of  the  complete  solid-propellant  rocket 
exhaust  are  not  necessary.  Small  mass  loadings  in  the  test  case 
are  therefore  required  to  simulate  the  larger  mass  loadings  of 
a  rocket  exhaust.  Small  mass  loadings  are  also  being  used  to 
allow  the  identification  of  individual  particle  behavior  and 
effects. 

The  test  section  materials  are  graphite,  graphite  phenolic, 
silica  phenolic,  and  tungsten.  Tungsten,  in  addition  to  being  of 
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interest  as  a  rocket  nozzle  material,  also  eliminates  the  effects 
of  surface  recession  and  chemical  reaction  in  assessing  particle- 
impact  effects.  The  test  section  configuration  is  a  3-inch-di ureter 
disk,  the  thickness  of  which  will  be  varied  between  1/4-  and  3/4- 
inch  to  allow  some  control  of  the  transient  temperature  response  of 
the  test  sections.  The  plasma-generator  exit  flow  impacting  on  the 
test  sections  will  be  circular  in  cross  section  with  a  diameter  of 
1/2-inch  or  slightly  greater. 

The  test  variables  are  the  test  section  material,  the  wall 
temperature,  the  plasma-generator  operating  conditions,  the  particle 
temperature,  the  particle  velocity,  the  particle  size,  the  particle 
mass  loading,  the  angle  of  incidence  of  the  particles  on  the  test 
section,  and  the  chemical  nature  of  the  exhaust  gas.  Some  of  these 
variables  are,  of  course,  interdependent  and  all  are  discussed 
below.  The  test  section  materials,  except  for  tungsten,  are  the 
same  as  those  used  in  the  chemical-erosion  program,  Graphitite  GX 
graphite  (Basic  Carbon  Corporation) ,  MX4500  graphite  phenolic 
(Fiberite  Corporation)  and,  MX2600  silica  phenolic  (Fiberite  Corpo¬ 
ration)  . 

The  wall  surface  temperature  is  important  in  terms  of  chemi¬ 
cal  effects,  both  between  particle  and  wall  and  test  gas  and  wall, 
and  can'  be  expected  to  influence  the  impact  behavior  of  the 
particle.  A  range  in  wall  temperatures  which  covers  those  repre¬ 
sentative  of  aft  closure,  throat,  and  expansion  cone  is  anticipated. 
Wall  temperature  control  will  be  accomplished  by  variation  of 
test  section  thickness  as  discussed  above,  changes  in  plasma- 
generator  operating  conditions,  and  variation  of  the  exposure  time 
of  the  test  section  to  the  particle-free  exit  stream  before  parti¬ 
cles  are  introduced. 

The  particle  temperature  is  tied  directly  to  the  plasma- 
generator  operating  conditions.  The  residence  time  of  the  particles 
in  the  plenum  chamber  will  be  sufficient  for  temperature  equili¬ 
bration  with  the  test  gas.  The  anticipated  particle  temperature 
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range  is  2000' K  to  3500'  K.  The  particle  velocity  is  dictated 
by  the  plasma-generator  operating  conditions  and  by  the  exit 
nozzle  size.  Particle  velocities  up  to  6000  ft/sec  are  antici¬ 
pated.  The  particle  mass  loading  will  be  varied  from  0  to  2  per¬ 
cent  of  the  total  exit  mass  flow,  but  most  tests  will  be  performed 
at  low  mass  loadings  to  allow  individual  particle  effects  to  be 
assessed.  The  angle  of  incidence  between  the  particles  and  the 
test  section  surface  will  be  varied  from  10°  to  40° .  As  dis¬ 
cussed  previously,  tests  will  be  performed  with  a  chemical 
reactive  test  gas,  Mixture  4,  and  an  inert  test  gas,  Mixture  5, 
both  of  which  closely  duplicate  the  thermodynamic  properties  of 
the  solid-propellant  rocket  combust j'  *.  products.  The  inert  gas 
allows  the  assessment  of  particle  effects  in  the  absence  of 
chemical  reaction,  whereas  the  chemically  reactive  gas  simulates 
both  the  chemical  and  thermodynamic  environment. 

3.3.3  Facility  modification  and  instrumentation 

The  general  facility  setup  in  preparation  for  the  particle- 
impact  -program  has,  for  the  most  part,  been  completed  in  that  the 
similar  activities  in  preparation  for  the  chemical-erosion  program 
are  applicable  here.  The  additional  facility  modifications  and 
the  instrumentation  techniques  are  discussed  below. 

A  new  plenum  extension  ring,  shown  in  Figure  35,  will  be 
used  for  injection  of  particles  in  the  plenum  chamber  of  the  arc- 
plasma  generator  downstream  of  the  arc  zone.  This  ring  also 
accommodates  the  secondary  gas  injection.  There  are  eight  particle 
injection  ports  uniformly  spaced  to  help  insure  proper  distribution 
of  particles  ir.  the  exit  flow.  The  particles  are  injected  in  a 
direction  counter  to  the  rotation  imparted  to  the  primary  gas  by 
the  arc  rotation.  This,  together  with  the  secondary  gas  injection, 
cancels  the  arc-induced  rotation  and  also  helps  to  insure  good 
mixing  and  a  uniform  particle  distribution. 


An  existing  particle  hopper  will  be  used  for  particle 
injection  andvf  low  control.  This  unit  is  capable  of  handling 
powders  with  a  size  distribution  for  which  the  average  particle 
size  is  as  low  as  5  microns.  The  hopper  has  been  checked  out 
with  the  12  micron  particles  for  the  range  of  mass  loadings  to 
be  used  and  found  to  perform  properly.  Operation  with  5  micron 
particles  is  possible  only  at  very  low  mass  loadings.  Particle 
metering  and  control  is  accomplished  by  a  variable  speed  vibrator. 
The  vibrator  causes  the  particles  to  fall  from  a  storage  canni- 
ster,  through  an  orifice,  and  into  a  carrier-gas  pickup  chamber. 
The  carrier  gas,  in  this  case  a  portion  of  the  secondary  injec¬ 
tion  gas,  delivers  the  particles  to  the  plenum  chamber  in  eight 
separate  feed  tubes  as  discussed  above. 

Due  to  the  nature  of  the  program,  a  number  of  specialized 
instrumentation  techniques  are  required  in  addition  to  the  stan¬ 
dard  data- acquisition  methods.  A  water-cooled  calorimeter  will 
be  used  in  place  of  the  test  section  to  determine  the  average 
cold-wall  heat  flux  over  the  central  3/4-inch-diameter  surface 
of  the  test  sections.  The  calorimeter  will  have  a  3-inch  guard 
ring  enclosing  its  3/4-inch  sensing  surface  to  duplicate  the  test 
section  geometry.  The  radiative  contribution  to  the  total  heat 
flux  will  be  assessed  with  a  HyCal  asymptotic  radiometer  with  a 
15°  view  angle  and  a  calibrated  output  to  500  Btu/ft2-sec.  Some 
of  the  test  sections  vill  contain  thermocouple  instrumentation 
to  assist  in  assessing  the  test  section  response  and  the  wall 
surface  temperature.  The  instrumentation  technique  will  be  the 
same  as  that  used  in  the  chemical-erosion  program. 

The  response  of  the  test  slab  to  the  simulated  rocket 
environment  will  be  photographed  with  a  high-speed  Traid  motion 
picture  camera  at  a  film  speed  of  1000  frames  per  second.  Real 
time  color  motion  pictures  and  post-test  still  photographs  will 
also  be  taken  to  further  evaluate  test  section  response.  In  some 
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of  the  tests,  additional  analyses  of  the  slab  response  and 
particle  behavior  will  be  performed  with  an  electron  microscope 
for  particle  identification,  sizing,  and  counting.  This  service 
will  be  performed  by  Allied  Space  Products,  Inc.,  of  Palo  Alto, 
under  Vidya  supervision. 

Particle-velocity  measurements,  both  before  and  after 
impact,  will  be  made  with  an  Abtronics  image-intensifier  camera 
as  discussed  in  Reference  2.  With  this  technique,  the  particle 
velocities  are  measured  by  a  streak  photography  method  in  which 
individual  particles  should  be  identifiable.  The  capability  of 
the  camera  system  in  performing  this  task  was  checked  out  further 
during  the  report  period.  Calculations  were  made  of  the  energy 
that  is  available  from  a  single  particle  for  energizing  the 
camera  system,  and  of  the  accuracies  obtainable.  The  calcula¬ 
tions  indicated  a  marginal  probability  of  the  camera  being  able 
to  "see"  the  particles  and  allowing  a  determination  of  particle 
velocity  from  the  resultant  photograph  with  an  error  no  grenfcer 
than  about  +5  percent.  Since  the  capabilities  do  appear  marginal, 
it  has  beer,  decided  to  perform  the  preliminary  tests  using  the 
camera  on  a  short-term  demonstration  basis.  In  the  preliminary 
tests,  only  a  one-dimensional  view  of  the  particle  stream  will 
be  taken  in  order  to  eliminate  a  complicated  split-image  optics 
system.  A  second  Traid  high-speed  motion  picture  camera  will 
also  view  the  same  field  as  the  image-iniensifier  camera  to 
obtain  a  continuous  record  of  the  particle  stream  both  before 
and  after  impact. 

An  indication  of  particle  temperatures  will  be  attempted 
with  the  two-color  pyrometer  discussed  previously.  Measurements 
will  be  made  both  before  and  after  impact  to  assess  the  energy 
lost  in  impaction  by  those  particles  that  are  reflected  from  the 
surface.  The  use  of  an  infrared  pyrometer  is  also  anticipated. 
The  pyrometer  measurements  will  be  made  on  the  stream  both  with 
and  without  particles  to  checkout  the  validity  of  the  pyrometer 
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indications.  Color-densitometer  analyses  of  photographs  of  the 
particle-laden  stream  as  discussed  in  Reference  2  will  also  be 
made  to  determine  particle  temperatures. 

Particle  sampling  for  determination  of  the  size  distribu¬ 
tion  in  the  high-temperature  exit  flow  of  the  plasma  generator 
still  remains  a  problem.  An  estimate  of  the  size  distribution 
of  particles  adhering  to  the  test  samples  will  be  made  as  part 
of  the  electron  microscope  examination  of  the  slab  surface. 

This,  of  course,  may  not  be  representative  of  the  free  stream 
distribution  due  to  the  possibility  of  a  selective  adherence 
of  the  particles  to  the  surface;  particles  in  one  size  range 
may,  for  the  most  part,  adhere  to  the  surface,  whereas  in  another 
size  range  they  may  tend  to  reflect.  0  .her  sampling  techniques 
are  under  investigation  and  include  sampling  by  collecting  the 
complete  exit  flow  over  a  period  of  time  and  "quick  action" 
sampling  uisng  a  nylon  disk  or  right-circular  cylinder.  In  the 
latter  case,  particles  are  collected  by  adherence  to  the  nylon 
pickup  which  has  been  softened  by  the  high-temperature  environ¬ 
ment  itself.  In  this  method  the  exposure  of  the  pickup  to  the 
stream  would  be  for  a  very  short  time  only,  1  second  or  less. 

This  technique  will  be  attempted  in  the  preliminary  test  program 
and  the  results  analyzed  to  determine  the  validity  of  the  distri¬ 
bution  obtained  from  it. 

3.3.4  Test  preparations 

Preparations  for  the  actual  performance  of  the  test  program 
were  begun  during  the  report  period.  Fortunately,  much  of  the 
setup  and  checkout  of  the  program  had  already  been  completed,  in 
that  much  of  that  required  for  the  particle-impact  program  is 
common  to  the  chemical-erosion  program.  The  progress  in  the 
additional  facility  and  instrumentation  preparations  are  discussed 
below. 
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The  new  plenum  extension  ring  was  designed  (Pig.  35)  and 
fabrication  begun.  The  particle-feed  unit  was  checked  out  and 
found  to  operate  satisfactorily  for  the  range  of  test  conditions 
to  be  used  in  the  program.  The  test  section  materials  for  the 
preliminary  program  are  on  hand.  The  excess  materials  from  the 
nozzle  fabrication  in  the  chemical-erosion  program  are  being  used. 
The  required  test  section  holder  and  actuating  system  are  a  part 
of  the  arc-plasma  generator  facility  and  required  no  modification 
to  accommodate  the  program. 

The  necessary  arrangements  for  use  of  the  Abtronics  image- 
intensifier  camera  on  a  demonstration  basis  were  made.  The  basic 
camera  will  be  used  with  "breadboard"  controls  to  allow  variation 
of  the  camera  operation  to  obtain  an  optimum  camera  system  for 
photographing  the  particles.  Overhaul  and  checkout  of  the  Traid 
high-speed  motion  picture  camera  recently  obtained  through  the 
government  surplus  list  for  use  on  this  program  was  begun.  The 
second  Traid  camera  is  ready  for  operation.  The  flat-face  heat- 
flux  calorimeter  is  a  facility  item  and  therefore  ready  for  use. 
The  guard  ring  around  it  required  to  duplicate  the  test  section 
geometry  will  be  fabricated  shortly.  The  HyCal  radiometer  has 
been  placed  on  order  for  delivery  by  the  end  of  the  first  month 
of  the  quarter. 

Preliminary  testing  for  checkout  of  the  complete  setup 
will  take  place  late  in  the  first  month  of  the  quarterly  period. 
The  preliminary  test  program  is  expected  to  be  completed  at  the 
end  of  the  first  week  of  the  second  month. 

The  future  effort  will  be  dictated  by  the  results  of  the 
preliminary  test  program.  If  no  unusual  problems  arise,  the 
particle-impact  program  is  expected  to  be  continued  to  completion. 
Otherwise  the  chemical-erosion  program  will  be  resumed  while 
appropriate  fixes  are  made. 
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TABLE  I.-  EXAMPLE  OF  OUTPUT  FROM  THE  ABLATING  SURFACE  THERMO¬ 
CHEMICAL  EQUILIBRIA  COMPUTER  PROGRAM. 


ABLATING  SURFACE  THERMO-CHEMICAL  EQUILIBRIA 


PRESSURE • ATM  2.72000E  00  SURF. TEMPeK  2745.48  ENTH.CAL/GM-4.41839E  01 

MOOT  P.G./CM  2.00000E-01  MOL  .WEIGHT  29.1972  HWCWM/CM1  -2.297S6E  02 

MOOT  CHAR/CM  A. 00000 E  00  SPEC.HEAT  0,31396  SURFACE  C2SI*  S 

-  -  -  -  MOLE  FRACTION  -  SPECIE 

0.0000000  CH  0.0000007  CHN  0.0000165  CHO  0.0000000  CH2 

0.0000000  CH3  0.0000000  CH4  0.0000000  CN  0.0221163  C02 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0143897  H  0.0000006  HN  0.0014957  HO 

0.0000000  HS1  0.0000000  HAS!  0.0000004  N  0.0001035  NO 

0.0000000  N02  0.0000000  NSI  0.0000000  N20  0.000 10^  0 

0.0000161  02  0.0000046  SI  0.0000000  S22  0.0019607  02SI 

0.3958704  CO  0.1256270  H2  0.0469780  H20  0.0853577  N2 

0.3059578  OS I 

PRESSUREtATM  2.72000E  00  SURF.TEMP»K  2718.20  ENTH»CAL/GM-6. 32286E  01 

MOOT  P.O./CM  2. 00000 E- 01  MOL .WEIGHT  30.6408  HW< 1 AM/CMI  -4.55246E  02 

MOOT  CHAR/CM  6.00000E  00  SPEC.HEAT  0.29580  SURFACE  02SI*  5 

•  -  -  -  MOLE  FRACTION  -  SPECIE 

0.0000000  CH  0.0000008  CHN  0.0000154  CHO  0.0000000  CH2 

0.0000000  CHS  0.0000000  CH4  0.0000000  CN  0.0169409  C02 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0117938  H  0.0000004  HN  0.0008443  KO 

0.0000000  HSI  0.0000000  H4S!  0.0000003  N  0.0000542  KO 

0.0000000  N02  0.0000000  NSI  0.0000000  N20  0.0000595  0 

0.0000064  02  0.0000064  SI  0.0000000  $12  0.0015461  02SI 

0.4249789  CO  0.1032327  H2  0.0272893  H2G  0.0647067  N2 

0.3485236  OS! 

PRESSUREtATM  2.72000E  00  SURF.TEMP.K  2644.91  ENTH,CAL/GM~9. 52692E  01 

MOOT  P.G./CM  2.00000E-01  MOL. WEIGH!  32.1360  HW(1*M/C*1  -1.06701E  03 

MOOT  CHAR/CM  1.00000E  01  SPEC.HEAT  0,27851  SURFACE  02SI*  5 

•  -  -  -  MOLE  FRACTIOM  -  SPECIE  -  - 

0.0000000  CH  0.0000013  CHN  Q.C000129  CHO  0.0000000  CH2 

0.0000000  CM3  0.0000000  CH4  0.0000000  CN  0.0090698  CG2 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0077538  H  0.0000002  HN  0.0002499  HO 

0.0000000  HSI  0.0000000  H4SI  0.0000001  N  0.0000142  NO 

0.0000000  N02  0.0000000  NS7  0.0000000  N20  0.0000154  0 

0.0000008  02  0.0000103  SI  0.0000000  SI2  0.0507956  02SI 

0.4575991  CO  0.0782426  H2  0.0100161  H20  0.0436389  N2 

0.3925796  OS I 


ABLATING  SURFACE  THERMQ-CmEM ICAl  EOUILIHKIA 


PRESSURE, ATM  2.720006  00  SURF*  T6MP»K  2191. 50  ENTH»CAL/GP  5.95612c  00 
MOOT  P.G./CM  4. 00000E-01  MOL. WEIGHT  16.7465  HW(5*M/CM1  8.39813E  00 

MOOT  CHAR/CM  10. 00000E-03  SPEC. HEAT  0.53985  SURFACE  G2SI*  6 

-  MOLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000019  CHN  0.0000076  CHC  0.0000000  CH2 

0.0000000  CH3  0.0000000  CH4  0.0000000  CN  0.0152161  C02 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0020477  H  0.0000000  HN  0.0000786  HO 

0.0000000  HSI  0.0000000  H4SI  0.0000000  N  0.0000024  NO 

0.0000000  N02  0.0000000  NS  I  0.0000000  N2U  0.0000003  0 

0.0000000  02  0.0000000  SI  0.0000000  Si 2  0.0000045  U2S1 

0.2923518  CO  0.4000892  H2  0.1079133  H20  0.1806630  N2 

0.0016238  OSI 

PRESSURE, ATM  2.72000E  00  SURF. f EMP, K  2324.41  ENTH,CAL/GM  7.81842E  01 
MOOT  P.G./CM  4.00000E-01  MOL. WEIGHT  16.9731  HWiWM/C*)  1.13367602 

MOOT  CHAR/CM  5.00000E-02  SPEC. HEAT  0.53626  SURFACE  02SI*  4 

-  MCLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000018  CHN  0.0000102  CHO  0.0000000  CH2 

0.0000000  CHi  0.0000000  CH4  0.0000000  CN  0.0143096  C02 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0041486  H  0.0000001  HN  0.0001940  HO 

0.0000000  HSI  0.0000000  H4SI  0.0000000  N  0.0000070  NO 

0.0000000  N02  0.0000000  NSI  0.0000000  N20  0.0000013  0 

0.0000001  02  0.0000000  SI  0.0000000  SI2  0.0000257  02S1 

0.2957008  CO  0.3931121  H2  0.1064358  H20  0.1780540  N2 

0.0079986  OSI 

PRESSURE, ATM  2.72000E  00  SURF. TEMP, K  2385.53  ENTH,CAL/GM  1.06596E  02 
MOOT  P.G./CM  4.00000E-01  MOL. WEIGHT  17.2575  HWIWM/CM)  1.59893b  02 

MOOT  CHAR/CM  1.00000E-01  SPEC. HEAT  0.52861  SURFACE  02SI*  4 

-  MOLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000018  CHN  0.0000116  CHO  0.0000000  CH2 

0.0000000  CH3  0.0000000  CH4  0.0000000  CN  0.0139588  C02 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2rt 

0.0000000  C4H  0.0055637  H  0.0000002  HN  0.0002804  HO 

0.0000000  HSI  0.0000000  H4SI  0.0000000  N  0.0000110  NO 

j  O.OOCOOOO  N02  0.0000000  NSI  0.0000000  N20  0.0000026  0 

0.0000003  02  0.0000000  SI  0.0000000  SI2  0.0000533  02SI 

:  0.2991946  CO  0.3861750  H2  0.1040248  H20  0.1750012  N2 

I  0.0157202  OSI 

{ 


=i 

I 


ABLATING  SURFACE  THERMO-CHEMICAL  EQUILIBRIA 


PRESSURE  *  ATM  2.T2000E  00  SURF. TEMP, K  24*7.66  ENTH,CAL/GM  1.27244E  02 
MOOT  P.G./CM  4.000C0E-01  MOL. WEIGHT  17.8068  HWIUM/CM)  2.03590E  02 

MOOT  CHAR/CM  2.00000E-01  SPEC. HE  AT  0.51295  SURFACE  C2SI*  4 

-  MOLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000018  CHN  0.0000152  CHU  0.0000000  CH2 

0.0000000  CH3  0.0000000  CH4  0.0000000  £N  0.0156321  C02 

9.0000000  C2H2  0.0000000  C5  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0073328  H  0.C000002  HN  0.0003932  HO 

0.0000000  HSI  0.0000000  H4SI  0.0000000  N  0.0000165  NU 

0.0000000  N02  0.0000000  NS  I  0.0000000  N20  0.0000050  0 

0.0000005  02  0.0000001  SI  0.0000000  SI2  0.0001078  02SI 

0.3056521  CO  0.3739429  H2  0.0992097  H20  0.1692827  N2 

0.0304091  0S1 

PRESSURE • ATM  2.72000E  00  SURF. TEMP 2508.25  ENTH.CAL/GM  1.32274E  02 
MOOT  P.G./CM  4.00000E-01  MOL  .WEIGHT  (8.8166  HW(1*M/CM)  2.38093E  02 

MOOT  CHAR/CM  4.00000E-01  SPEC. HE AT  0.48522  SURFACE  02SI*  4 

-  -  MOLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000018  CHN  0.0000149  CHO  0.0000000  CH2 

0.0000000  CH3  0.0000000  CH4  0.0000000  CN  0.0(32976  C02 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0093411  H  0.9000003  HN  0.0005193  HO 

0.0000000  HSI  0.0000000  H4SI  0.0000001  N  0.0000231  NO 

0.0000000  N02  0.0000000  NS I  0.0000000  N20  0.0000088  l) 

0.0000009  02  0.0000004  SI  0.0000000  SI2  0.0002064  02SI 

0.3173353  CO  0.3527691  H2  0.0903566  H20  0.1590020  N2 

0.0571220  0S1 

PRESSURE. ATM  2.72000E  00  SURF. TEMP»R  2541.11  ENTH. CAL/GM  1.24593E  02 
MOOT  P.G.SCM  4.00000E-01  MOL.  WEIGH  I  19.7176  HWd+M/CMI  2.49185E  02 

MOOT  CHAR/CM  6.00000E-01  SPEC. HEAT  0.46239  SURFACE  Q2SI*  4 

-  -  •  -  MOLE  FRACTION  -  SPECIE  -  -  -  - 

0.0000000  CH  0.0000019  CHN  0.0000159  CHO  0.0900000  CH2 

0.0000000  CH3  0.0000000  CH4  0.0000000  CN  0.0130509  C02 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  9.0104806  H  0.0000004  HN  0.0005828  HO 

0.0000000  HSI  0.0000000  H4SI  0.0000001  N  0.0000268  NO 

0.0000000  N02  0.0000000  NSI  0.0000000  N20  0.0000116  0 

0.0000012  02  0.0000006  SI  0.0000000  SI2  0.0002896  02SI 

0.3277504  CO  0.3344700  H2  0.0825563  H20  0.1499513  N2 

0.0800094  0S1 


ABLATING  SURFACE  THERMO-CHEMICAL  EQUILIBRIA 


PRESSURE, ATM  2.72000E  00  SURF • TEMP, K  2576.90  ENTH*  C AL /GM  1.02127E  02 

MOOT  P.G./CM  4.00000E-01  MOL. WEIGHT  21.2527  HW(1*M/CM)  2.45105E  02 

MOOT  CHAR/CM  I.OOOOOE  00  SPEC. HEAT  0.42758  SURFACE  C2SI*  4 

-  MULE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000019  CHN  0.0000171  CHO  0.0000000  CH2 

0.0000000  CH3  0.0000000  CH4  0.0000000  CN  0.0125714  C02 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0116111  H  0.0000004  HN  0.0006225  HO 

0.0000000  HSI  0.0000000  H4SI  0.0000001  N  0.0000295  NU 

0.0000000  N02  0.0000000  NSI  0.0000000  N20  0.0000150  0 

0.0000014  02  0.0000012  SI  0.0000000  SI2  0.0004144  02S1 

0.3456076  CO  0.3038944  H2  0.0695342  H20  0.1346853  N2 

0.1209921  CS1 

PRESSURE, ATM  2. 72000E  00  SUKF. TEMP,K  2607.67  FNTH,CAL/GM  5.19700E  01 
MOOT  P.G./CM  4.00000E-01  MOL. WEIGHT  24.0231  HWll+M/CM)  1.76698E  02 

MOOT  CHAR/CM  2.00000E  00  SPEC. HEAT  0.37556  SURFACE  02SI*  5 

-  MOLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000021  CHN  0.0000178  CHO  0.0000000  CH2 

0.0000000  CH3  0.0000000  CH4  0.0000000  CN  0.0112015  C02 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

2000000  C4H  0.0119348  H  0.0000004  HN  0.0005447  HO 

G. 0000000  HSI  0.0000000  H4SI  0.0000001  N  0.000026 2  NU 

0.0000000  N02  0.0000000  NSI  0.0000000  N2U  0.0000163  0 

0.0000013  02  0.0000028  SI  0.0000000  S12  0.0005592  02SI 

0.3784511  CO  0.2495631  H2  0.0470340  H20  0.1074637  N2 

0.1931805  OSI 

PRESSURE, ATM  2.72000E  00  SURF. TEMP, K  2601.00  ENTH.C AL/GM- 1 . 12858E  01 

MOOT  P.G./CM  4.00000E-01  MOL. WEIGHT  27.2046  HWCl+M/CM)  -6.09432E  01 

MOOT  CHAR/CM  4.00000E  00  SPEC.HE AT  0.32855  SURFACE  02SI*  5 

-  MOLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000025  CHN  0.0000169  CHO  0.0000000  CH2 

0.0000000  CH3  0.0000000  CH4  0.0000000  CN  0.0081892  C02 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0100664  H  0.0000003  HN  0.0003013  HO 

0.0000000  HSI  0.0000000  H4SI  0.0000001  N  0.0000140  NO 

0.0000000  N02  0.0000000  NSI  0.0000000  N2U  0.0000102  0 

0.0000005  02  0.0000059  SI  0.0000000  SI2  0.0005244  02SI 

0.4177502  CO  0.1874090  H2  0.0233317  H2U  0.0766250  N2 

0.2757523  OSI 


ABLATING  SURFACE  THERMO-CHEMICAL  EQUILIBRIA 


PRESSURE, ATH  2.72000E  00  SURF. TEMP*K  2563.09  EN f H» C AL /GM-5. 315876  01 

MOO?  P.G./CM  4.000006-01  MOL. WEIGHT  28.9663  HW(1*M/CM)  -3.93374E  02 

MOOT  CHAR/CM  6.00000E  00  SPEC. HE AT  0.306SS  SURFACE  C2Sl»  6 

-  MCLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000035  CHN  0.0000151  CHO  O.OOOOOOC  CH2 

0.0000000  CH3  0.0000000  CH4  0.0000000  CN  0.0053497  C02 

0.0000000  C2H2  0.0000000  C5  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0077620  H  0.0000002  HN  0.0001360  HO 

0.0000000  HSI  0.0000000  H4SI  0.0000001  N  0.0000060  NO 

0.0000000  *02  0.0000000  NS  I  0.0000000  N20  0.0000044  0 

0.0000001  02  0.0000093  SI  0.0000000  SI2  0.0003615  U2SI 

0.4409863  CO  0.1523724  H2  0.0115624  H20  0.0595792  N 2 

0.3218520  OS  I 

PRESSURE, ATM  2.72GOOE  00  SURF. TEMP,K  2336.22  ENTH.CAL/GM- 1 .4984  IE  02 

MOOT  P.G./CM  4.00000E-0I  MOL. WEIGHT  30.9482  HWCHM/CM1  -1.70819E  05 

MOOT  C  1AR/CM  1.00000E  01  SPEC .HE A!  0.28413  SURFACE  02SI*  7 

-  -  MCLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000301  CHN  0.0000090  CHO  0.0000000  CH2 

0.0000000  CH3  0.0000000  CH4  0.0000000  CN  0.0005646  C02 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0023992  H  0.0000000  HN  0.0000029  HO 

0.0000000  HSI  0.0000000  H4SI  0.0000000  N  0.0000001  NO 

0.0000000  N02  0.0000000  NS I  0.0000000  N20  0.0000000  0 

0.0000000  02  0.0000309  SI  0.0000001  SI2  0.0000304  02S1 

0.4683762  CO  0.1143621  H2  0.0007770  H20  0.0412801  N2 

0.3721371  OSI 

PRESSURE, ATM  2.72000E  00  SURF. TEMP, K  2072.65  6NTH, C AL/GM  1.75800E  02 

MOOT  P.G./CM  6.00000E-01  MOL. WEIGH!  15.3372  HWd+M/CM)  2.83038E  02 

MOOT  CHAR/CM 10.00000E-03  SPEC. HE A!  0.56095  SURFACE  02SI*  10 

-  -  -  -  MOLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000087  CHN  0.0000971  CHO  0.0000000  CH2 

0.0000000  CH3  0.0000002  CH6  0.0000000  CN  0.0047066  C02 

0.0000000  C2H2  ■  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0010941  H  0.0000000  HN  0.0000086  HU 

0.0000000  HSI  0.0000000  H4SI  0.0000000  N  0.0000002  NO 

0.0000000  NC2  0.0000000  NS  I  0.0000000  N20  0.0000000  0 

0.0000000  02  0.0000000  SI  0.0000000  SI2  0.0000008  02SI 

0.3382361  CO  0.4779035  H2  0.0318263  H2Q  0.1449022  N2 

0.0013053  OSI 


ABLATING  SURFACE  THERMO-CHEM I CAL  EQUILIBRIA 


PRESSURE, ATM  2.720006  00  SURF. TEMP, K  2189.83  ENTH, CAL /GM  2.34S9G6  02 
MOOT  P.G. /CM  6.000006-01  MOL. WEIGHT  15.5411  HWIUM/CM)  3. 87074b  02 

MOOT  CHAR/CM  5.00000E-02  SPEC. HEAT  0.55756  SURFACE  C2SI»  4 

-  MOLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000086  CHN  0.0000095  CHQ  0.0000000  CH2 

0.0000000  CH5  0.0000001  CH4  0.0000000  CN  0.0043670  C02 

0.0000000  C2HZ  0.0000000  Ci  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0022032  H  0.0000000  HN  0.0000208  HO 

0.0000000  HSI  0.0000000  H4SI  0.0000000  N  0.0000005  NO 

0.0000000  N02  0.0000000  NS  I  0.0000000  N20  0.0000001  0 

0.0000000  02  0.0000000  SI  0.0000000  SI2  0.0000044  02 S 1 

0.3402471  CO  0.4720524  H2  0.0313649  H20  0.1432691  N2 

0.0064523  OSI 

PRESSURE, ATM  2.72000E  00  SURF. TEMP, K  2243.20  ENTH, CAL/GM  2.53158E  02 
MOUT  P.G. /CM  6. OOOOOE-0 1  MOL. WEIGHT  15.7944  HWIUM/CM1  4.30368E  02 

MOOT  CHAR/CM  1.000006-01  SPEC.HE Af  0.55025  SURFACE  02SI»  4 

-  MCLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000086  CHN  0.0000108  CHU  0.0000000  CH2 

0.0000000  CH5  0.0000001  CH4  0.0000000  CN  0.0042059  C02 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0029462  H  0.0000001  HN  0.0000298  HO 

0.0000000  HSI  0.0000000  H4SI  0.0000000  N  0.00C0008  NO 

0.0000000  N02  0.0000000  NSI  0.0000000  N2Q  0.0000001  0 

0.0000000  02  0.0000001  SI  0.0000000  S 12  0.0000091  C2S 1 

0.3425529  CO  0.4656185  H2  0.0305663  H20  0.1413220  N2 

0.0127288  OSI 

PRESSURE, ATM  2.72000E  00  SURF. TEMP, K  2296.86  ENTH, CAL/GM  2.594706  02 
MOOT  P.G. /CM  6. OOOOOE-O 1  MOL. WEIGHT  16.2852  HW(1+M/C»M  4.67046E  02 

MOOT  CHAR/CM  2.00000E-01  SPEC. HEAT  0.53513  SURFACE  02SI*  4 

-  MOLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000087  CHN  0.0000121  CHO  0.0000000  CH2 

0.0000000  CH3  0.0000001  CH4  0.0000000  CN  0.0040133  C02 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0038690  H  0.0000001  HN  0.0000410  HO 

0.0000000  HSI  0.0000000  H4S I  0.0000000  N  0.0000012  NO 

0.0000000  N02  0.0000000  NSI  0.0000000  N20  0.0000002  0 

0.0000000  02  0.0000002  SI  0.0000000  SI2  0.00001A2  02S1 

0.3469550  CO  0.4537429  H2  0.0289304  H20  0.1376177  N2 

0.0247898  OSI 


ABLATING  SURFACE  THERMO-CHEN ICAL  EQUILIBRIA 


PRESSURE* ATM  2.T2000E  00  SURF.TEMP.K  2347.94  ENTH.CAL/GM  2.43597E  02 

MOOT  P.G./CM  6.00000E-01  MOL. WEIGHT  17.1994  HWIHM/CM)  4.87195E  02 

MOOT  CHAR/CM  4.00000E-01  SPEC. HEAT  0.507B4  SURFACE  02SI*  4 

-  MOLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000090  CHN  0.0000135  CHO  0.0000000  CH2 

0.0000000  CH3  0.0000000  CH4  O.OUOOOOO  CN  0.0037464  C02 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

0.0000000  C4H  0.0048968  H  0.0000001  HN  0.0000521  HG 

0.0000000  HSI  0.0000000  H4S I  0.0000000  N  0.0000016  NO 

0.0000000  N02  0.0000000  NS  I  0.0000000  N20  0.0000004  0 

0.0000000  02  0.0000005  SI  0.0000000  SI2  0.0000342  02S1 

0.3551371  CO  0.4323386  H2  0.0258350  H20  0.130B079  N2 

0.0471267  OSI 

PRESSURE* ATM  2.72000E  00  SURF. TEMP, K  2374.31  ENTH.CAL/GM  2.198176  02 

MOOT  P.G./CM  6.00000E-01  MOL. WEIGHT  18.0308  HWIUM/CM)  4.835986  02 

MOOT  CHAR/CM  6.00000E-01  SPEC. HE AT  0.48488  SURFACE  02S1*  4 

-  -  -  -  MOLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000094  CHN  0.0000143  CHU  0.0000000  CH2 

0.0000000  CH3  0.0000000  CH4  0.0000000  CN  0.0035197  C02 

0.0000000  C2H 2  0.0000000  C3  0.0000000  C3H  0.0000000  C2H 

O.OOOOOOG  C4H  0.0054509  H  0.0000001  HN  0.0000560  HO 

0.0000000  HSI  0.0000000  H4SI  0.0000000  N  0,0000018  NO 

0.0000000  N02  0.0000000  NS  I  0.0000000  N20  0.0000005  0 

0.0000000  02  0.0000009  SI  0.0000000  SI2  0.0000468  02SI 

0.3624033  CO  0.4132301  H2  0.0230298  H20  0.1246642  N2 

0.0673722  OSI 

PRESSURE  * ATM  2.72000E  00  SURF.TEMP«K  2399.56  ENTH*CAL/GM  1.724796  02 

MOOT  P.G./CM  6.00000E-01  MOL. WEIGHT  19.4846  HW(1*M/CM)  4.484446  02 

MOOT  CHAR/CM  1.00000E  00  SPEC. HEAT  0.44897  SURFACE  02SI*  4 

-  -  -  -  MOLE  FRACTION  -  SPECIE 

0.0000000  CH  0.0000103  CHN  0.0000149  CHO  0.0000000  CH2 

0.0000000  CH3  0.0000000  CH4  0.0000000  CN  0.0030927  C02 

0.0000000  C2H2  0. ) 000000  C3  0.0000900  C3H  0.0000000  C2H 

0.0000000  C4H  0.0059051  H  0.0000001  HN  0.0000538  HO 

0.0900000  HSI  0.0000000  H4SI  0.0000000  N  0.0000017  NO 

0.0000000  N02  0.0000000  NS!  0.0000000  N20  0.0000005  0 

0.0000000  02  0.0000018  SI  0.0000000  SI2  0.0000627  02SI 

0.3757439  CO  0.3802465  H2  0.0181962  H20  0.1139891  N2 

0. 1026807  OSI 


ABLATING  SURFACE  THERMO-CHEMICAL  EQUILIBRIA 


PRESSURE, ATM  2.72000E  00  SURF. TEMP, K  2401.18  ENTH,CAl/GM  7.88366E  01 

MOOT  P.G./CM  6.00000E-01  MOL. WEIGHT  22.2368  HW( l+M/CMI  2.83812E  02 

»  MOOT  CHAR/CM  2.00000E  00  SPEC.HSA1  0.39309  SURFACE  Q2SI*  5 

-  MOLE  FRACTION  -  SPECIE  - 

0.0000000  CH  0.0000149  CHN  0.000C146  CHO  0.0000000  CH2 

0.0000000  CH3  0.0000000  CH4  0.0000000  CN  0.0020346  C02 

0.0000000  C2H2  0.0000000  C3  0.0000000  C3H  0.0000000  C2M 

0.0000000  C4H  0.0054481  H  0.0000001  HN  0.0000307  HO 

0.0000000  HSI  0.0000000  H4SI  0.0000000  N  0.0000010  NO 

0.0000000  N02  0.0000000  NSI  0.0000000  N20  0.0000003  0 

I  0.0000000  U2  0.0000048  SI  0.0000000  SI2  0.0000639  02SI 

0.4010091  CO  0.3187139  H2  0.0094089  H2U  0.0939510  N2 

0.1693039  OSI 

PRESSURE, ATM  2.72000E  00  SURF. TEMP,*  2206.32  ENTH.CAL/GM-8. 1 3879E  01 

MOOT  P.G./CM  6.00000E-01  MOL. WEIGHT  25.6498  HWIHM/CM)  -4.55772E  02 

4  MOOT  CHAR/CM  4.00000E  00  SPEC. HEAT  0.33842  SURFACE  SIC  9 

-  MOLE  FRAC1I0N  -  SPECIE  - 

0.0000000  CH  0.0006017  CHN  0.0000091  CHU  0.0000000  CH2 

0.0000012  CH3  0.0000041  CH4  0.0000000  CN  0.0000498  C02 

0.0000089  C2H2  0.0000000  C3  0.0000000  C3H  0.0000001  C2H 

0.0000000  C4H  0.0017391  H  0.0000000  HN  0.0000001  HQ 

0.0000001  HSI  0.0000001  H4SI  0.0000000  N  0.0000000  NO 

0.0000000  N02  0.0000000  NSI  C.0000000  N20  0.0000000  0 

I  0.0000000  02  0.0000939  SI  0.0000032  SI2  0.0000016  02S1 

0.4328691  CO  0.2440134  H2  0.0001467  H20  0.0693722  N2 

0.2510864  OSI 

PRESSURE, ATM  2.72000E  00  SURF. TEMP»K  2595.75  ENTH,CAL/GM  3.53855E  01 

MOOT  P.G./CM  6.Q0000E-01  MOL. WEIGHT  27.5572  HWCUM/CM)  2.68930E02 
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Humber  in  parenthesis  indicates  initial  nozzle  throat  diameter  in  inches 
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Number  in  parenthesis  indicates  initial  nozzle  throat  diameter  in  inches. 
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Number  in  parenthesis  ii  jHcates  initial  nozzle  throat  diameter  in  inches 


Figure  1.-  Program  flow  chart. 


erosion  calculation 
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Figure  3.-  Theoretical  chevical  erosion  of  graphite  in  an 
H,  C,  M,  0,  C\  exhaust. 


Figure  4.-  Typical  X-ray  photograph  for  determining 

thermocouple  locations. 
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.  (a)  0 . 3-inch-diaswter  throat  nozzle,  high- 

teeperature  condition. 


Figure  5.*-  Pressure  distribution,  Mixture  4. 


Pressure  ratio,  p/p 


nssaJa 


I 


(c)  0. 5-inch-dia»eter  throat  nozzle,  high- 
tamper  attire  condition. 


Figure  5.-  Continued 
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(d)  0.4-inch-diameter  throat  nozzle,  low 
temperature  condition. 


Pigure  5.-  Concluded. 
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(a)  0.3-inch-diameter  throat  nozzle,  high- 
temperature  condition. 


Figure  6.-  Prasaure  distribution.  Mixture  5. 
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(b)  0.4-inch-diameter  throat  nozzle,  high- 
temperature  condition. 
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Figure  6.-  Continued. 


Pressure  ratio,  p/p 
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(c)  0. 5-inch-diameter  throat  nozzle,  high- 
temperature  condition. 


Figure  6.-  Continued. 
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(a)  0.3-inch-diameter  throat  nozzle,  high- 
temperature  condition. 

gure  7.-  Heat -transfer  coefficient  axial  distribution,  Mixture  4. 
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(b)  0. 4-inch-aiametar  throat  nozzle,  high- 
temperature  condition. 

Figure  7.-  Continued. 


Heat -transfer  coefficient,  lb/ft  -sec 


(c)  0 . 5 -inch -diameter  throat  nozzle,  high- 
temperature  condition. 


Figure  Concluded. 


HeaL-transfer  coefficient.  Ik-'f-'  sec 


1.0 


0.6 


0.6 


0.4 


0.  2 


0 


(a)  0.3-inch-diameter  throat  nozzle,  high- 
temperature  condition. 


Figure  8.-  Heat-transfer  coefficient  axial  distribution.  Mixture  5. 
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(b)  0. 4-inch-diameter  throat  nozzle,  high- 
temperature  condition. 


Figure  8.-  Continued. 
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(c)  0, 5-inch-diameter  throat  nozzle,  high 
temperature  condition. 


Figure  8.-  Concluded. 
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Figure  10.-  Cut-away  of  fired  nozzle,  Graphitite  GX  graphite,  Mixture 
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Figure  10.-  Continued. 
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figure  10,  -  Concluded 
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Figure  11.-  Cut-away  of  fired  nozzle,  ATJ  graphite, 
Kixture  4.  Nozzle  B43,  nominal  conditions  of 
T0  -  3500°  K  and  fp)  -  300  paia,  9£  -  30.2  sec, 

Test  Ho.  1000.  &m° 
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Figure  12.-  Cut-away  of  fired  nozzle,  MX4500  ophite  phenolic,  Mixture 
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Figure  12.-  Concluded. 
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Figure  13.-  Cut-away  of  fired  nozzle,  MX2600  silica  phenolic,  Mixture 
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Figure  14.-  Cut-away  of  fired  nozzle,  Graphitite  GX  graphite,  Mixture 
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Figure  16.-  Put- away  of  fired  nozzle,  MX4500  ophite  r>,  »ol.ir  iixture 
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Figura  16.-  Concludad. 
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(c)  Nozzle  E25,  nominal  conditions  of 

Tq  *  2000°  K  and  (pQ)  -  300  psia, 

0f  *  29,7  sec>  Te8t  No.  1031. 
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(a)  Nozzle  E17,  Test  No.  939. 

Figure  18. -  Surface  recession  and  surface  recession  rate, 
Graphitite  GX  graphite,  Mixture  4. 
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(b)  Vossle  110,  Test  Vo.  940. 
Figure  18.  -  Continued. 
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(c)  Nozzle  E19,  Teat  No.  997. 
Figure  18.-  Continued. 
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Figure  19.-  Surface  recession  and  surfaco-recession  rate, 
ATJ  graphite,  Mixture  4,  Nozzle  E43,  Test  No.  1000. 
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Figure  20.-  Surface  recession  and  surface-recession  rate,  MX4500 
graphite  phenolic,  Mixture  4,  nozzle  B38,  Text  Mo.  1001. 
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Figure  21.- Surface  recession  and  surface-recession  rate.  MX2f>00  silica  phenolic. 

Mixture  4.  Nozzle  E22.  Tes*  No.  941. 


MX2600  silica  phenolic 


Figure  22.-  Surface  recession  and  surface-recession  rate,  MX2600  silica  phenolic, 

Mixture  5,  nozzle  E24,  Test  Ho.  960. 
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Figure  25.-  Comparison  of  surface  recession  with  and  without  swirl. 
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(a)  Nozzle  E10,  Test  Mo.  940. 

Figure  27.-  Yeaqperature  histories,  Graphitite  GX  graphite.  Mixture 


Figure  2  7.-  Concluded. 


Temperature  histories,  MX4500  graphite  phenolic 
Mixture  4,  Nozzle  E22,  x’est  Mo.  941. 


Figure  29. -  Temperature  histories,  MX2600  silica  phenolic. 
Mixture  4,  nossle  122,  Test  Vo.  941. 


Qraphitita  A?  graphite 


HQ  «  3100  Btu/lb 

T  «  3600°  K 
o 

pQ  «  177  psia 

Dc  *  0,399  inch 
Thermocouple  locations 

rx  «  0.353  inch 
va  »  0.753  inch 
r  -  1.159  inch 


Te^c  time,  sec 

(b)  Mossle  *32,  Test  *o.  1004 


Figure  30.-  Concluded 
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Figure  31.-  Teaperature  histories,  MX4500  graphite  phenolic. 
Mixture  5,  Nozzle  E32,  Test  Vo.  1004. 


Temperature 


Figure  32.-  Temperature  histories.  MX2600  silica  phenolic. 
Mixture  5.  nozzle  E32.  Test  Ho.  1004. 


Figure  33.-  Char  recession  and  char  depth,  MX2600  silica  phenolic 

Mixture  4,  nozzle  E22,  Test  Mo.  941. 
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Pigu_e  34.-  Char  recession  and  char  depth.  MX2600  silica  phenolic. 

nozzle  E24,  Test  No.  960. 


